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ABSTRACT 

We employ the Chandra Multiwavelength Project (ChaMP) and the Sloan Digital Sky Survey 
(SDSS) to study the fraction of X-ray-active galaxies in the field to z = 0.7. We utilize spectro- 
scopic redshifts from SDSS and ChaMP, as well as photometric redshifts from several SDSS catalogs, 
to compile a parent sample of more than 100,000 SDSS galaxies and nearly 1,600 Chandra X-ray 
detections. Detailed ChaMP volume completeness maps allow us to investigate the local fraction of 
active galactic nuclei (AGN), defined as those objects having broad-band X-ray luminosities Ljs:(0.5- 
8 keV) > 10^^ ergs~^, as a function of absolute optical magnitude. X-ray luminosity, redshift, mass, 
and host color/morphological type. In five independent samples complete in redshift and i-band ab- 
solute magnitude, we determine the field AGN fraction to be between 0.16 ± 0.06% (for z < 0.125 
and -18 > Mi > -20) and 3.80 ± 0.92% (for z < 0.7 and Mi < -23). We find excellent agreement 
between our ChaMP/SDSS field AGN fraction and the Chandra cluster AGN fraction, for samples 
restricted to similar redshift and absolute magnitude ranges: 1.19 ± 0.11% of ChaMP/SDSS field 
galaxies with 0.05 < z < 0.31 and absolute i?-band magnitude more luminous than Mr < —20 are 
AGN. Our results are also broadly consistent with measures of the field AGN fraction in narrow, 
deep fields, though differences in the optical selection criteria, redshift coverage, and possible cosmic 
variance between fields introduce larger uncertainties in these comparisons. 
Subject headings: galaxies: active; galaxies: nuclei; surveys; X-rays: galaxies 



1. INTRODUCTION 

A fundamental constraint on all theories modeling the 
interplay of supermassive black hole (SMBH) accretion 
and galaxy evolution should be the fraction of galaxies in 
the local universe that host actively accreting nuclei. Ac- 
cretion onto a SMBH is the predominant source of energy 
produced by active galactic nuclei (AGN) and QSOs, so 
a measure of accretion activity (at a variety of redshifts 
and luminosities) provides a valuable cons traint on the 
black hole mass function at the present day ()Soltanlll982l : 
Eecs 1984: lMarcom et al.l[200l . 

In both active and inactive galaxies, observations show 
a remarkably tight correlation between the properties 
of galactic spheroids and the mas ses of their central 
black holes, the Mbh — cr relation (iFerrarese fc MerrittI 
[2000t IGebhardt et"all [2OOOI: IGiiltekin et all I2009D . sug^ 
gesting that SMBHs play an important role in the for- 
mation of galactic bulges. Simulations indicate that the 
black hole and the bulge may co-evolve, each growing 
as a result of repeate d mer ger events (e.g., JHcrnau isti 
IT989I: iPTMatteo et all [20051: iHopkins et all 120051 120081 1 
In this scenario, the merger causes cold gas to fall 
to the center, triggering starbursts, and fueling rapid 
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black hole growth. Even while accretion inflows increase 
the black hole's mass, outflows transport energy from 
the active nucleus t o the surrounding host galaxy (e.g., 
i Granato et"all l200l. creating a mechanism by which the 
galaxy and it's central black hole can "feedback" as they 
evolve. In extreme cases, these outflows can carry enough 
energy to alter the galaxy's evolution by ejecting gas and 
truncating star formation. 

Several well-studied phenomena support this so- 
called feedback paradigm, which links star formation 
and SMBH growth to galaxy mergers. First, the 
space density of low-luminosity active galactic nuclei 
(LLAGN) pea ks at lower host m a ss than that of lumi- 
nous quasars i Barger e t all 120051 : iHasinger et al.l 120051 : 
iScannapieco et al. 2005^^ indicating that AGN may 
undergo "cos mic downsizing" as do galaxy spheroids 
(jCowie et al.l [l996iV Next, the local mass density of 
SMBHs corresponds closely to the lu minosity density 
produced by quasars at high redshift ()Yu fc Tremaind 
120021 : [Hopkins et all 120061 ). Hence, the present-day pop- 
ulation of AGN and quiescent black holes may represent 
the fossil record of their brighter, high-redshift counter- 
parts (the QSO luminosity function peaks at redshifts 
z ~ 2), formed at a time when galaxy mergers were com- 
mon. This suggests that all massi ve galaxies may have 
undergone an AGN phase (e.g., IMarconi et al.l [200^ 
IShankar et"an [2004). 

Though they clearly influence the evolution of galax- 
ies and the conditions in the early universe, fundamental 
questions about the formation and evolution of accreting 
black holes remain. To be robust, the merger and feed- 
back paradigm must reproduce the fraction of galaxies 
hosting AGN, the ratio of obscured to unobscured AGN 
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necessary to match the cosmic X-ray background, the 
space density of SMBHs, and AGN clustering. 

X-ray emission is the most rehable primary signature 
of AGN activity because X-rays offer the most complete 
and efficient marker of active accretion close to the black 
hole (fO-fOO gravitational radii). Combined with other 
wavelengths, X-ray observations sample a wide variety 
of AGNs and improve our census of AGN demographics 
(|Brandt fc Hasingeill2005l and references therein). 

Galaxies that do not harbor a powerful accreting black 
hole may also produce X-rays due to hot interstellar gas, 
supernova remnants, a nd low- and high-mass X-r ay bi- 
nary populations ( Fabbianol 120061 : iKim et al.ll2006D : and 
even in "normal" galaxies a LLAGN may persist (see lHol 
120081 for a review). In this work, we consider the X- 
ray-active fraction of galaxies out to z = 0.7, in a field 
covering more than 20 deg^, at a wide range of broad- 
band X-ray luminosities (10'*° < [0.5-8 keV; ergs"*] 
< lO**^), including those characteristic of normal and 
starforming galaxies. Since the many X-ray-producing 
mechanisms described above can muddy the waters in 
lower-X-ray-luminosity systems, we focus the bulk of 
our analysis on sources with X-ray luminosities Lx(0.5- 
8 keV) > 10^^ ergs^*, characteristic of accreting SMBHs. 

Measurements of the AGN fraction as a function 
of environment have been made using both optical 
spectroscopic samples a nd op tical /radio-selected sam- 
ples. iKaufF mann et al.l (|2004l ) use emission -line galax- 
ies in the Sloan Digital Sky Survey (SDSS; lYork et all 
120001 ) to measure the fraction as a function of envi- 
ronment and find that the fraction of very luminous 
AGN decreases with increasing environmental density 
(see also |Popcsso_&__Biviano 200^; while less luminous 
SPSS AGN' show no such trend (|Miller et al.l [200l 
iKauffmann et al.ll2004l ). Studies of the radio-loud SDSS 
AGN fraction demonstrate that the fraction of galaxies 
hosting radio-loud AGN increases towards richer environ- 
ments (except a t the highest emission-line luminosities, 
iBest et al.ll2005[) . 

As in field studies. X-ray identification of AGN activity 
in clusters has proved to be even more efficient than op- 
tical or radio techniques, particularly for low-luminosity 
AGN. A significant body of work on the cluster AGN 
fraction has been enabled by a combination of Chandra 
X-ray observations and spectroscopic follo w-up for a se- 
ries of 32 c lusters with 0.05 < 2 < 1.3 (jMartini et al.i 
[200l[200l[2?M[200^ These authors find larger AGN 
fractions in clusters than those found in optical emission 
line studies. One of the aims of the present work is to 
make a detailed comparison between the field and cluster 
X-ray AGN fractions based on a combination of Chandra 
imaging and ground-based optical phot ometry. 

NA SA's Chandra X-ray Observatory (jWeisskopf et al.l 
|2002| ) has generated a rich archive of X-ray sources, de- 
tected with unprecedented resolution and sensitivity. By 
cross-correlating the Chandra archive with the SDSS — 
an extensive optical survey in five filters — it is possible 
to identify X-ray counterparts to individual optical ob- 
jects, confirming their classification as quasars or active 
galactic nuclei. The Chandra M ultiwavelength Project 
(ChaMP: [Green et al.|[200ll2009[ ) has carefully analyzed 
323 Chandra fields (about 30 square degrees) that overlap 
the SDSS and characterized all optical/X-ray matches. 



We employ these X-ray and optical surveys to determine 
the fraction of actively accreting SMBHs (or AGN) out 
to redshift of ~ 0.7. We also explore the AGN fraction as 
a function of galaxy luminosity, redshift, mass, and host 
galaxy morphology. Such measures of the AGN fraction 
will prove useful for constraining the fueling, lifetimes, 
and growth of central SMBHs. 

In ^J2] we describe the Extended ChaMP {Chandra 
Cycles 1-6) and SDSS surveys; in ^J3] we detail the 
ChaMP/SDSS sample developed for this study, including 
specifics on the spectroscopic and photometric redshift 
catalogs we employ. Section 2] outlines our methods for 
determining k-corrections, absolute magnitudes, masses. 
X-ray and optical luminosities, and our definition of five 
independent samples complete in redshift and i-band ab- 
solute magnitude. In Sj5] we describe our determination 
of the X-ray-active and AGN fractions in these volume- 
limited optical samples, investigate possible evolution of 
the fraction with redshift, and trends with absolute mag- 
nitude and mass, in addition to providing an assessment 
of how photometric redshift errors and other selection ef- 
fects might impact our measure of the AGN fraction. We 
discuss the AGN fraction as a function of galaxy color 
and compare our results to others in the field, groups, 
and clusters in ^ and wrap up with our conclusions in 

m 

We adopt an Hq = 70 kms"*Mpc~*, = 0.7, and 
= 0.3 cosmology throughout. 

2. THE X-RAY AND OPTICAL CATALOGS 

2.1. The Extended Chandra Multiwavelength Project 

The Chandra Multiwavelength Project is a wide-area 
serendipitous X-ray survey based on archival X-ray im- 
ages of the |6| > 20 deg sk y observed with the A XAF 
CCD Spectrometer (ACIS; IWeisskopf et"aI1 120021) on- 
board ChandrE- The full 130-fie ld Cycle 1-2 X-ray 
catalog is public (jKim et al.ll2007al ). and the most com- 
prehensive X-ray number counts (log iV-log S) to date 
have been produced thanks to 6600 sources an d massive 
source- retrieval simulations (jKim et al.ll2007bl ). 

We have recently expanded our X-ray analysis to cover 
a total of 392 fi elds through Chandra Cycle 6 (see deta iled 
descriptions in lCovev et al.l [20081 : iGreen et al.ll2009| ). to 
improve statistics and encompass a wider range of source 
types — we refer to this expansion as the Extended 
ChaMP. The new list of Chandra pointings avoids over- 
lapping X-ray observations by eliminating the observa- 
tion with the shor ter exposure time. As described in 
IGreen et all (pOOl . we also avoid Chandra fields with 
large (>3 ') extended sources in either optical or X- 
rays (e.g., nearby galaxies MlOl, NGC 4725, NGC 4457, 
or clusters of galaxies MKW8, or Abell 1240). Spuri- 
ous X-ray sources (due to, e.g., hot pixels, bad bias, 
bright source readout streaks) have been fiagged and re- 
moved as described in IKim "etaI[([2007aD . For the expan- 
sion, we select only Chandra fields that overlap SDSS 
Data Release 5 imaging; of the 392 ChaMP fields (ob- 
servation IDs; hereafter "obsids"), 323 fall within the 
SDSS DR5 footprint The Exten ded ChaMP covers 
-30 deg'^ at the brightest fiuxe^ ([Green et al.l [2009). 

^ ChaMP results and data are available on-line: 
http: / /hea- www. harvard.edu/CHAMP 

" A list of the Chandra fields included in this study (ob- 
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In addition to SDSS imaging, we have undertaken deep 
(r ^ 25) NOAO /MO SAIC optical imagining (described 
in lGreen et al.|[2004D . work currently being extended to 
67 fields (W.A. Barkhouse et al. 2010, in preparation). 

We have incorporated spectroscopic information, 
redshifts and line diagnostics, for the entire opti- 
cal catalog fr om: (1) existing ChaMP spectroscopy 
(|Green et al.l [200l: (2) the SPSS DR7 spectro- 
scopic catalog ()Abazaiian et al.l I2009D : (3) the SDSS 
DR7 Max Planck Institute for Astrophysics/ John's 
Hopkins University (MPA/JHU) value-added galaxy 
catalog ([Brinchmann et aTl l2004fF^ : and (4) cross- 
correlat ion with the NA SA Extragalactic Database 
(NED) . iGreen et al.l ([20091 ) and iConstantin eTall ([2009f ) 
give details on matching the optical spectroscopy with 
the Extended ChaMP X-ray source catalog. 

The Extended ChaMP also includes a comprehensive 
set of sensitivity maps for ACIS imaging, implemented 
in the xskycover table. This allows (1) recognition of 
imaged-but-undetected objects, (2) counts limits for 50% 
and 90% detection completeness, (3) corresponding flux 
upper limits at any sky position, as well as (4) flux sen- 
sitivity versus sky coverage for any subset of obsids. To 
generate the ta ble we use the wavde tect detection algo- 
rithm in CIAO ([Freeman et al.ll2002l ) to generate thresh- 
old maps at each wavdetect kernel scale (1, 2, 4, 8, 16, 
and 32 pixels). The threshold maps, computed from the 
local background intensity, determine the magnitude of 
the source counts necessary for a detection at each pixel 
with a detection threshold of P = 10~^ (corresponding 
to 1 false source per 10^ pixels). Thus, when an object is 
not detected at a given location, the threshold map value 
serves as an upper limit to the source counts, as normal- 
ized via th e detail ed simulated source retrieval results of 
iKim et all ([2007al ). The xskycover table covers the full 
Extended ChaMP with sky pixels, each 10' x 10', whose 
boundaries are chosen to match a regular commensurate 
grid across the sky. The final sensitivity in any given sky 
pixel is interpolated from the two threshold maps com- 
puted at scales that bracket the size of the point spread 
function (PSF) at th at location (for ad ditional details, 
see the Appendix to IGreen et al.l l2009| l . The accuracy 
of the xskycover counts limits for 50% and 90% detec- 
tion completeness were verified bv lAldcroft eFall ([2001 
using the Chandra Deep Field South (CDF-S). 

2.2. The Sloan Digital Sky Survey 

The Sloan Digital Sky Survey is an extensive 
photometric and spectroscopic optical survey that 
covers n early o ne-qu a rter of the northern sky 
fStouahto n et al.l 120021: lAbazai ian et al. 2003, 2004, 
12005; .Adelman-McCarthv et al.ll200&i2007i .i2008). With 
a dedicated 2.5m telescope the SDSS accomplishes 
uniform photometry in five filters (w, r, i, z) to magni- 
tudes as faint as r < 22.5, with astrometric uncertainty 

sid, R.A., Dec, exposure time [ksec], number of ChaMP/SDSS 
matches, and total number of optical spectra) can be found at 
http: / /hea- www. cfa.harvard.edu/CHAMP/ 

^ The brightest ChaMP fluxes are ~ fx{0.5 - S.OkeV) ~ 5 X 
10~^'' ergcm~^s~^, where the flux limit represents the number 
of counts detectable in 90% of simulation trials, converted to flux 
assuming a power-law F = 1.7 at z = and the Galactic Nh 

appropriate for each obsich 

Publicly available at |http://www.mpa-garching.mpg.de/SDSS/| 



similar to Chand ra^ s (i.e., better than 0.1" for r < 20.5; 
iPier et al.l I2003D . The SDSS data pipeline analyzes 
object morphology a nd provides r eliable star-galax y sep- 
arati on to r ^ 21.5 ([Lupton et a l. 2002; Scranton et al.l 
1200^ . During its main survey phase (through Data 
Release 7 [DR7]), the SDSS has obtained photometric 
measurements for over 350 million unique objects. The 
SDSS also employs a pair of dual multiobject fiber-fed 
spectrographs, and has obtained spectra for ^900,000 
galaxies, ^120,000 quasars, and -450,000 stars 0. 

For the present study we query SDSS Data Re- 
lease 5 (DR5: lAdelman-McCarthv et ar]i2007t ). i.e. imag- 
ing and spectroscopy taken through July 1, 2005, 
which includes photometry for 217 million objects over 
8000 deg^ and about 1 million spectra of galaxies, 
quasars, and stars chosen from 5713 deg'^ of the imag- 
ing data. By selecting SDSS data within 20' of each 
Chandra aimpoint (or within 28' where the expan- 
sion yielded additional high-confidence matches), we 
create a ChaMP/SDSS catalog of objects within the 
323 Chandra obsids described above. A photomet- 
ric recalibration effort for the SDSS Data Release 6 
(DR6), resulted in in iprovements to the SDSS pho- 
tometric zero points dAdelman-McCarthv et al.l 120081 : 
iPadmanabhan et al.ll2008D . We thus cross-correlate our 
DR5 ChaMP/SDSS catalog with the DR6 to utilize these 
"UberCal" magnitudes. As mentioned above, we also 
cross-correlate our ChaMP/SDSS catalog with the new 
SDSS DR7 spectroscopic database to obtain the largest, 
most up-to-date spectroscopic sample. 

3. SAMPLE SELECTION 

To create a high-quality optical galaxy sample with 
an X-ray flux (or flux limit) for each object, we ap- 
ply a series of X-ray and optical selection criteria to 
the ChaMP/SDSS sample. Our sample-selection logic 
is summarized in Fig. [1] and described in detail below. 

3.1. X-ray and Optical Quality Cuts 

As described in ^2.1\ the Extended ChaMP contains 
392 unique Chandra fields, 323 of which fall within the 
SDSS DR5 footprint. The catalog of SDSS objects re- 
covered within this overlap region contains 980,214 in- 
dividual SDSS objects; we refer to this as our "Parent" 
sample. 

To insure the deepest and most uniform X-ray cover- 
age, we exclude area from the ACIS-S S4 chip (CCD 8) 
which sustained damage early in the Chandra mission 
and has high background and streaking as a result. We 
accept only area with off-axis angles (OAA) < 12' — 
at larger off-axis angles the PSF broadens considerably 
and becomes distorted ('Fei gelson et al.l [20021 : IKim et all 
|2007b), which leads to large uncertainties in centroid- 
ing and source counts. We also exclude area near ACIS 
chip edges and in chip gaps where the exposure times 
are significantly shorter. To enforce this criteria, we re- 
quire that the effective exposure time at each xskycover 
sky pixel be greater than 80% of the maximum exposure 
time for the obsid. These X-ray quality cuts reduce our 
X-ray sky coverage to ~26 deg . 

See |http://www.sdss.or g' for additional details and public data 
access. 
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Parent sample 
(980,214) 

SDSS objects from DR5 within 
ChaMP sl<ycoverage area 



X-ray quality cuts: 



1. NotonACIS-S4 

2. Off-axis angle < 12 amin 

3. Not on chip edge/gap 
(effective exposure time/max exposure time > 0.8) 



NO 



X-ray-to-Optlcal 
Matching: 


YES (483,940) 


1. High-quality, visually inspected, matched sources 

2. Not cosmic-ray X-ray afterglow (nor affected by one) 


Optical quality cuts: 


YES (480,529) 



NO 



Reject 

(51% rejected) 



Reject 

(0.7% rejected) 



1 . SDSS primary source 

2. SDSS star or galaxy 

3. Not saturated 

4. 14.5 < psfMagJ < 21.0 (extinction corrected) 

5. No SDSS photometry error codes (ugriz > 0) 

6. Bad u-amp runs cut (see Richards et al. 2002) 

7. High Proper Motion cut (see Munn et al. 2004, 2008) 



Redshift Measurement Selection 
and Redshift Cut: 



NO 



Reject 

(74% rejected) 



YES (123,509) 



YES 



Spectroscopy? 



NO 



ChaMP 
or, if not, 
SDSS 



YES 



0.0025 < z < 0.995? 



Galaxy photoz 



Spatially Resolved? 



QSO photoz 



YES 



NO 



QSO Candidate? 



Photoz & 
Color cuts 



YES 



NO 



General photoz 



NO 

Reject 
(15% rejected) 



YES (105,550) 



X-ray-detected? 




Parent X-ray-observed Sample 
(105,550) 



YES 



Parent X-ray-detected Sample 
(1,593) 



Fig. 1. — Decision tree showing the sample selection criteria for clean samples of SDSS objects within the ChaMP sky coverage area for 
which Chandra X-ray flux limits can be assigned, and for Chandra X-ray detections. We refer to these as the "Parent X-ray-observed" and 
"Parent X-ray-detected" samples, respectively. 
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In a number of the Chandra fields the SDSS DR5 data 
does not cover the full obsid. To crudely quantify the 
optical coverage of our X-ray fields, we perform visual 
inspection of the 323 ChaMP/SDSS overlap obsids and 
assign a coverage fraction to each. An estimate of the 
sky area based on these rough coverage fractions is ap- 
proximately 21 deg^ and is a lower limit to our actual 
joint ChaMP-SDSS coverage. The SDSS catalog associ- 
ated with this sky coverage area contains 483,940 objects 
(^49% of the Parent sample). 

For the X-ray-detected subset of the galaxy sample, we 
apply additional quality constraints. All X-ray-detected 
objects in the ChaMP survey area are visually inspected, 
wherein we overplot X-ray centroids and their associated 
positions on both the X-ray and optical images to iden- 
tify poorly matched, multiply matched, or photometri- 
cally contaminated objects. We accept only objects with 
the highest-confidence, uncontaminated matches, i.e. a 
single optical counterpart with an X-ray-to-optical posi- 
tion offset no greater than 2" and/or less than the 95% 
X-ray position error. Match statistics for the X-ray de- 
tections in the Extended ChaM P c atalog are describe d 
in detail in iCovev et all (|2008[ ) and iGreen et al.l (|2009f ). 
We also eliminate X-ray sources that most likely result 
from (or have possibly been affected by) cosmic-ray after- 
glows. The sources rejected by this set of ChaMP quality 
cuts decreases our sample by less than 1%. 

An accurate AGN fraction requires a clean sample of 
both X-ray-detections and non-detections. We thus also 
remove photometric contaminants from our optical sam- 
ple using the "flags" included in the SDSS photomet- 
ric catalog^. We include only unsaturated, PRIMARY 
(mode=l) objects that are classified as either a star 
(type= 3) or a galaxy (type= 6) and whose 5-band pho- 
tometry is without error codes (u, g,r,i, z > 0). To avoid 
poor photometry both at the bright and faint extremes of 
the survey, and to insure uniform coverage at the faintest 
fluxes, we select objects with i-band psf-magnitudes be- 
tween 14.5 and 21.0. We also eliminate objects impacted 
by bad w-amplifier runs (Richards ct al. 2002). 

As described in §2.2) the SDSS pipeline provides reli- 
able star-galaxy separation for objects down to r ~ 21.5. 
This distinction is a morphological one, and thus cannot 
correctly characterize active galaxies and quasars that 
appear to be point sources without consulting optical 
color information. Faint objects at high redshift are also 
likely to be misclassified by morphological diagnostics. 
Eliminating all SDSS "stars" would thus bias our overall 
fraction. As a result, we do not remove point sources 
a priori from our sample. Instead, we eliminate stars 
based on their proper motions, positions in SDSS color- 
color and color-magnitude diagrams, and spectroscopic 
or photometric redshifts. (See ^ 33. 21 below for a descrip- 
ti on of our color and redsh ift cuts.) 

iMunn et al.l (|2004 120081) have published a proper mo- 
tion catalog that combines USNO-B and SDSS observa- 
tions. We use the Munn catalog to assign proper motions 
to our catalog objects where there is a single USNO-B 
match within 1" (match= 1), the rms fit residuals (sigRa, 
sigDec) are less than 350 mas in each coordinate, and 
there is at least one USNO-B detection and one SDSS 

See http://www.sdss.org/dr7/products/cataIogs/flags.html 
for a detailed description ot the yUSB image processing Hags. 



0.20 



0.15 



0.10 



0.05 



0.00 




Fig. 2. — Fractional histogram of the best redshift for the "Full" 
sample (105,550 galaxies; black solid line) and the spectroscopic 
sub-sample (2668 objects from ChaMP and SDSS DR7 spec- 
troscopy; gray dashed line). The Full sample contains redshifts to 
^ ~ 1, where the photometric redshifts become less reliable, while 
the spectroscopic sub-sample shows a strong peak at low-redshift, 
0.05 < bestz < 0.1. The peak in the spectroscopic sub-sample 
results largely from SDSS spectroscopic selection criteria. 

detection per source (nfit > 2). To eliminate high proper 
motion objects (i.e. stars), we require that the proper 
motion in at least one coordinate be larger than 3 a, 
where a is the proper motion uncertainty in that coordi- 
nate (e.g., |PAfii;yi| > 3 X (JjiA), and that the total proper 
motion be larger than 11 masyr^^. 

The combined X-ray and optical quality cuts described 
in this section result in a catalog of 123,509 SDSS objects. 

3.2. Redshift and Color Selection 

To assign a redshift to each object in our Parent sam- 
ple we exploit both spectroscopic and photometric red- 
shift catalogs. We select a bestz, the highest quality 
redshift available, by ranking the redshift catalogs as 
follow s: (1) the C haMP spectroscopic catalog (398 ob- 
jects; IGreen et all [20041: (2) the S DSS DR 7 spectro- 
scopic catalog (2,270 objects; Abaza iian et al.l 2009): (3) 
the SDSS neural network (NN) Photoz2 catalog fo r re- 
solved objects (68,908 objects: lOvaizu et al.l [20081) : (4) 
the DR6 Quasar Photometric Redshift Catalog (281 ob- 
jects; Richards ct al. 2009); and (5) for unresolved point 
sources with no other redshift determination, the DR6 
version of the SDSS Photoz catalog, derived from the 
template-fitting metho d with repaired inte rpolated tem- 
plates (33,693 objects: [Csabai et al.ll2003[ ). Properties 
of the two SDSS photometric redshift catalogs, Photoz 
and Photoz2, are summarized in Table [T] We choose 
the DR6 version of the SDSS Photoz catalog (hence- 
forth, the "point source template" or "template" photo- 
z) because it contains photometric redshift estimates for 
point sources that are not included in the other cat- 
alogs; e.g., the newer DR7 catalog is limited to ob- 
jects classified morphologically by the SDSS pipeline as 
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Fig. 3. — Spectroscopic redshift versus photometric redshift for 
2668 objects in our galaxy sample with spectroscopy (398 from 
ChaMP and 2270 from SDSS DR7). The symbol type represents 
the photometric redshift c atalog used: t he SDSS DR6 resolved 
NN photo-2 galaxy catalog (Ovaizu et al.''2008, black points), the 
SDSS DR6 QSO catalog (Richards ct al. 2009, blue filled circles 
and blue arrows), the SDSS DR6 point source template photo- 
z catalog l lCsabai et al.|[2003 . green stars). The most pro mine nt 
photo-2 failures at redshifts above 2 ~ 1 are for QSOs (see i|5.4l l. 
(A color version of this figure is available in the online journal.) 



galaxies (|Abazaiian et al.l l2009| ) . Based on our selec- 
tion scheme, only objects with stellar morphology are 
assigned point source template photo- 2;'s. We eliminate 
all spectroscopically-confirmed stars from our sample. 

Figure [2] shows normalized histograms of the redshifts 
for our "Full" sample (solid black line) and for the spec- 
troscopic sub-sample (gray dashed line). For those ob- 
jects with both photometric and spectroscopic redshifts, 
97% have < 0.1 (Fig. ^. The photo-z algorithms 

are the least reliable for quasars (Fig. [3l filled circles and 
arrows). We discuss the impact of these potentially large 
photometric redshift errors on our determination of the 
AGN fraction in 

To quote a single error estimate on our bestz value, 
we utilize and/or convert the redshift errors in each cat- 
alog into a 1 (J Gaussian error. ChaMP spectroscopic 
redshifts and 1 a redshift errors are estimated via radial 
velocities, using the cross-correlation an d emission line 
fitting techniques in the IRAF rvsao task (iKurtz fc MinkI 
119981: iGreen et al.ll2004l) . Redshifts and their errors for 
the SDSS DR7 spectroscopic sample are estimated using 
emissi on-line or cross- correlation techniqueJ*^. 

The [ Richards etlll §009) QSO photometric redshift 
catalog uses a Bayesian approach to estimate a most 
likely redshift {zphot), as well as a high and low red- 
shift {zphothi and zphotlo) which represent an approx- 
imate range for zphot, and a probability that zphot is 

A discussion of SDSS redshift determinations is available at 
|http://www. sdss.org/d r7/algorithms / rcdshift_type.html 



Catalog 


Photoz 


Photoz2 


Number 


33,693 


68,908 


Source type 


unresolved 


resolved 


Method 


template fitting 


neural network 


Abbreviation 


template 


NN 


Reference 


Csabai et al. (2003) 


Ovaizu et al. f2008) 









TABLE 1 
SDSS Photometric Redshifts 



within this range. This photo metric redshift determ i- 
nation is described in detail by iWeinstein et al.l ()2004l) . 
We estimate 1 a errors (i.e. 68% probability intervals) by 
approximating the zphot errors as symmetric and Gaus- 
sian. We assume that the zphothi and zphotlo values 
bracket the zphotprob range for a Gaussian distribution 
and derive a "1 cr" error using Gauss' formula (i.e. using 
the error function, erf). We also utilize the full redshift 
probability distribution function for the QSO subset in 
the Monte Carlo simulations described in ^5.41 

Both the SDSS resolved NN photo-z and point source 
template photo-z catalogs include Gaussian 1 a errors, 
though the methods used to derive them d i ffer. In the 
resolved NN photo-z analvsis. lOvaizu et al.l (|2008D use a 
nearest neighbors estimator, with a spectroscopic train- 
ing set, to determine their best photometric redshift and 
the associated error. We select their "CC2" values, which 
employ both magnitudes and concentration indices in 
the redshift determination. The point source template 
photo-z catalog uses a template-fitting algorithm with 
repaired interpolated templates to determine the photo- 
metric redshift a nd associated 68% u ncertainty, assuming 
Gaussian errors (ICsabai et al.l[2nnl . 

For the sample that has been assigned a point source 
template photo-z, we make additional color and magni- 
tude cuts to reduce our contamination from stars. We 
accept all objects for which {u — g) < 0.6 (Fig. 2^, cyan 
contours), as these are likely to be extragalactic objects 
associated with a star forming or otherwise active galax- 
ies (see, for example, iRichards et al.l[2002 :). For objects 
with (u — g) > 0.6, we remove objects from two regions 
in color-color and color-magnitude space that are popu- 
lated primarily by stars (solid black lines. Figs. \^ and 
[4]i): (1) the first elimination region is defined by SDSS 
psfMagJ < 18 and (r — i) < 0.2; (2) the second is an area 
along the main sequence stellar locus in the (r — i) vs. 
(i — z) diagram (parametrized bv lCovev et al.ll2007l black 
dashed line. Fig. HJ;). For the second cut, we approxi- 
mate the main sequence stellar locus as a straight line, 
described by the formula (r — i) = 1.8 x (i — z) +0.05. We 
remove objects with (r — i) > 0.75 in a diagonal region 
within ± 0.3 of the stellar locus center line (solid black 
lines. Fig. |4t). We do not cut blueward of (r — i) = 0.75, 
because our spectroscopic sample (open purple circles), 
which consists of galaxies and AGN only, indicates that 
this would begin to remove galaxies from our sample. 
Despite our efforts to eliminate stellar point sources, our 
sample is still likely to suffer from considerable stellar 
contamination. We discuss and quantify the impact of 
this contamination on our AGN fraction in 

Once a bestz has been assigned and the most likely 
stellar contaminants have been removed, we restrict our 
attention to those objects with 0.0025 < bestz < 0.995. 
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Fig. 4. — Color-color and color-magnitude diagrams for the Full sample after the cuts described in i|3.2l and outlined in Fig. \T\ have been 
applied to remove galactic (stellar) contamina nts, without el iminating extragalactic objects. Gray contours/dots correspond to objects 
assigned resolved N N photo-2's from t he Ovaizu et al.l II2008I ) catalog. Red and cyan contours represent objects that have point source 
template photo-2's IjCsabai et al.ll2b03l 1 a nd (u — q) color s abov e and below 0.6 (panel a, solid black line), respectively. Dark blue stars 
indicate objects with QSO photo-z's from [Richards et al.l l|2009l). Objects with spectroscopy are indicated with open purple circles. The 
solid black lines in panels c and d indicate regions from which red ((u— q) > 0.6) poin t sources with template photo-z's have been eliminated. 
The black dashed lines in panels a-c show the SDSS stellar locus from lCovev et all l [2007 ). 



This redshift cut avoids cases where the photometric er- 
rors are large and/or the photometric redshift algorithms 
cannot find a suitable fit. These redshift, color, and spec- 
troscopic cuts remove another ^15% of the sample, leav- 
ing 105,550 objects in our optical galaxy catalog — the 
"Parent X-ray-observed" sample. Of these, 1,593 are X- 
ray detected — these form our "Parent X-ray-detected" 
sample. (As a result of this selection scheme, the ob- 
jects in the X-ray-detected sample are a subset of those 
in the X-ray-observed sample.) These objects are shown 
in the color-color and color-magnitude diagrams of Fig. 



131 ChaMP/SDSS spectroscopic redshifts (open purple 
circles); QSO photometric redshifts (dark blue stars); re- 
solved NN photo-z galaxies (gray contours); the point 
source template photo-z population with {u — g) < 0.6 
(cyan contours); and the point source template photo-z 
population with (u — g) > 0.6, which also meet the color 
cuts described above (red contours). 



DERIVED PROPERTIES AND VOLUME-LIMITED 
SAMPLES 
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4.1. Galaxy Properties 

We calculate optical k-correc tions and absolute magn i- 
tudes using kcorrect v4_l_4 (jBlanton fc Rowei^l2007[ ). 
via the relation 

mR = Mq + DM{z) + Kqb{z) - 5 log /i, (1) 

where rriR is the apparent magnitude, Mq is the absolute 
magnitude, DM{z) = 25 + 5 log [DL/{h-^ Mpc)] is the 
distance modulus calculated from the l uminosity distanc e 
Dl, and Kqr{z) is the k-correction (jHogg et al.l 120021 ) . 
In particular, we use sdssJscorrect and input our 
6este, extinction-corrected SDSS model-magnitudes, and 
model-magnitude errors. Since the kcorrect algorithm 
assumes a Hubble constant = 100 h kms~^Mpc~^ 
with /i = 1, we also convert to h = 0.7. The result- 
ing i-band k-corrections (to rest-frame z = 0) and abso- 
lute magnitudes, M^, are shown as a function of bestz in 
Fig. [5j We identify spectroscopic absorption line galax- 
ies (red stars) and narrow emission line galaxies (cyan 
stars) based on their MPA/JHU line classifications. The 
spectroscopic sub-sample appears at the bright end of the 
absolute magnitude distribution, as expected from SDSS 
main galaxy sample target selection algorithm which 
requires an r-band Petrosian magnitude brighter than 
17.77. In the top panel of Fig. [5] the absorption line 
galaxies follow a tighter trend in k-correction than their 
emission line counterparts. This trend is also expected, 
since emission line strengths vary widely in star-forming 
galaxies and AGN. We indicate objects with QSO pho- 
tometric redshifts (blue filled circles) and note that in 
many cases their k-corrections (and resulting Mi) lie far 
from the main galaxy locus. We test the impact of this 
bias on our fractions in ii5.4l and discuss its source in ijS.SI 

In the bottom panel of Fig. [5] we color-code the points 
by the redshift catalog from which they originate (red 
and cyan symbols as in the top panel for spectroscopic 
redshifts, blue filled circles for QSO photo-z's, light gray 
points for the resolved NN photo-z catalog, and green 
points for the point source template photo-z catalog). 
There are several noticeable structures at the bright end 
of the Mi vs. hestz relation. The most prominent of 
these results from the Ca H-K 4000 A break passing be- 
twee n the SDSS g and r filter s (see, for example. Fig. 
1 of IPadmanabhan et al.l [20071 ). We estimate the edges 
of this "gap" at 5200 A and 5700 A and mark the red- 
shift corresponding to observed-frame 4000 A as black 
vertical bars. In this rendering of the data, we see that 
most of the structure comes from the point source tem- 
plate photo-z objects, which are, by definition, red point 
sources {{u—g) > 0.6) in the SDSS. That the point source 
template photo-z catalog does not perform well for these 
objects is not a surprise since its primary spectroscopic 
training sets are the SDSS main galaxy sample and the 
luminous red galaxy sample. 

Using the ChaMP's X-ray sensitivity maps (see de- 
scription of the xsky cover table in fJHJ), we assign 
a Chandra broad-band (0.5-8.0 keV) X-ray fiux limit 
to each of the galaxies in our X-ray-observed sam- 
ple. We also assign a broad-band X-ray flux to each 
detection. With these X-ray fluxes (and flux lim- 
its) , w e calcu late th e X-ray-to-optical flux ratio, a la 
iHornschemeier et al.l ((2000.) . using the SDSS r-band ap- 
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Fig. 5. — {Top) K-correction as a function of bestz for the Full 
galaxy sample (gray points). Objects with spectroscopic classifi- 
cations from the DR7 MPA/JHU catalog are marked (red stars 
for absorption line galaxies, cyan stars for narrow emission line 
galaxies), and objects with QSO photometric redshifts (blue filled 
circles) are also overplotted for comparison. (Bottom) Absolute 
magnitude as a function of bestz, determined after the k-correction 
from the top panel has been applied. The light gray points indicate 
objects assigned a resolved NN photo-z and the green points are 
those with a point source template photo-z (i.e., the sum of the 
red and cyan contours in Fig. [4J; the remaining symbols are as in 
the top panel. The prominent structures at bestz ~0.3 and 0.43 
result from the 4000 A Ca H-K break passing between the SDSS 
g and r filters (see, e.g., [Padmanabhan et al. 200'3). We estimate 
the edges of this filter gap at 5200 A and 5700 A and mark the 
redshift corresponding to observed-frame 4000 A as black vertical 
lines. 

parent model-magnitudes: 

log (/.//.') - log (0.5 -2kcV) + 5.57 +(r72.5) , (2) 

where the conversion constant is determined assuming 
f{vo) = 3.631 X 10~^° ergcm"^ s~^ H z~^ and using 
the S DSS r'-band effective wavelengths (|Fukugita et al.l 
[1991). We also use the bestz, together with ChaMP X-ray 
fluxes and flux limits, to calculate luminosity distancef^ 
{Dl), broad-band X-ray luminosities (Lx) for detec- 



tions, and X-ray luminos ity limits (Lx 



90% confl- 

dence in the broad band, lAldcroft et al.l [2008) for non- 
detections. In Fig. [6| we show log Lx (dark gray stars) 
and log Lx.iim (light gray points) in units of ergs~^. On 
the left, the X-ray luminosities (limits) are shown as a 
function of bestz — the X-ray-detections span the locus 
of optical galaxies uniformly, indicating that our optical 
selection has not introduced obvious selection biases in 

^'^ Luminosity and comoving distances are calculated using the 
Python version of the Wright Cosmology Calculator (Wright 2006.). 
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bestz log [0.5-3.0 keV, erg s"'] 



Fig. 6. — (Left) Log of the broad-band (0.5-8.0 keV) X-ray luminosity, log Lx, as a function of bestz lov the X-ray-observed (light gray 
points) and X-ray-detected samples (dark gray stars). (Right) Absolute j-band magnitude vs. log Lx (gray symbols as in the left panel). 
The solid line marks log Lx = 42, the limiting X-ray luminosity used for our AGN fraction calculations; the dashed line indicates log 
Lx = 42.5, the limit used for AGN fraction calculations in the zCOSMOS survey (Silverman et al. 2009a). The black points (upper limits) 
and large, black open circles (detections) mark objects that reside off the main locus in the Mi vs. host mass diagram (see Fig. 1141 and 
i|5.5p . Representative 1 a error bars are plotted in the lower right corner of both panels. These are the median uncertainties for the 
X-ray-observed sample, derived from the 100-iteration "full" Monte Carlo simulation described in i|5.4l — the error bars thus derive only 
from uncertainties in the redshifts. 



the X-ray. On the right, we show Mi vs. log Lx (or log 
Lx,iim) and mark two luminosity limits, log Lx — 42 
and log Lx — 42.5, which we employ later in calculat- 
ing our field AGN fraction and comparing our results to 
those of deep-field surveys. 

Assuming an optical continuum power-law slope of 
a = —0.5 {fn oc z^", where v is the emission frequency), 
we derive rest-frame monochromatic optical luminosities 
at 2500 k{l • ; ergs"! Hz"!) using the SDSS g-band 

^ 2500A\ ° / o n 

dereddened magnitude. In theory the u-band magni- 
tude would be more appropriate for objects with hestz < 
0.53, because the filter's central wavelength is closer to 
(1 + z) X 2500A, but the large errors on u-band magni- 
tudes mitigate this improvement. We derive the 2 keV lu- 
minosity (Z2kcv; ergs^^Hz^^) from the broad-band flux 
using PIMMS and assuming an unabsorbed power-law 
with spectral slope F — 1.7. The slope of a hypothetical 
power law from 2500 A to 2 keV is commonly character- 
ized by aox, where Uox = 0.3838 x log .V^2kov); we 
plot aox as a function of I ^ in Fig. [71 

4.2. Volume- Limited Samples 

Using these derived properties, we define five sam- 
ples complete in bestz and Mi (referred to as volume- 
limited samples hereafter) to calculate our X-ray-active 
and AGN fractions. We note that this makes ours 
an optically-selected. X-ray-detected (or limited) study. 
The five volume-limited samples shown in Fig. [5] are 
defined as: (1) z < 0.125 and -18 > M, > -20; (2) 
z < 0.275 and -20 > M, > -21; (3) z < 0.4 and 
-21 > M, > -22; (4) z < 0.55 and -22 > M, > -23); 



I 'I 




0.5 - . - 

I , ■ ■ I I I I 

26 27 28 29 30 31 

log l,s„ (erg s"' Hz"'] 

Fig. 7. — aox vs. optical (2500 A) luminosity for the Full galaxy 
sample. X-ray-observed galaxies are shown as light gray points. 
X-ray-detections are marked with large symbols according to their 
broad-band X-ray-luminosity: log Lx > 42 (black diamonds) and 
log Lx < 42 (dark gray stars). The ChaMP QSO best-fit regres- 
sion line (aox = (0.061 ± 0.009)log(« ? ) - (0.319 ± 0.258)) from 

^^^^^^^^^^^^^^^^^^ 

I Green et al.l II2009I ) is shown as a black solid line with errors. Rep- 
resentative 1 (7 error bars are as described in Fig. |6] the median 
error bar for the aox is smaller than the width of the horizontal 

line for the log(i • ) error bar. 

''^ 2500A' 

and (5) z < 0.7 and M, < -23. (The low-redshift limit 
for all samples is z = 0.0025.) Their constituents are 
detailed in Tabled 

We determine each galaxy's stellar mass via the 
kcorrect tool, which generates the coefficients of a fit 
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TABLE 2 

Volume-Limited Samples 



s 


bestz 
max 


Mi 






Full 




Spec 


mm 


max 


xdet ^ xlim 




1 


0.125 


-18 


-20 


29 


5028 


2 


182 


2 


0.275 


-20 


-21 


88 


10949 


7 


342 


3 


0.400 


-21 


-22 


239 


17984 


37 


566 


4 


0.550 


-22 


-23 


416 


22095 


68 


567 


5 


0.700 


-23 




323 


13628 


57 


464 



Note. — Properties of the five volume-limited samples for the Full 
galaxy catalog (described in ^4.2| and shown in Fig. [Sj and the Spec- 
troscopic sub-sample. The first column gives the sample number; the 
second gives the maximum rcdshift [bestz; all samples have a mini- 
mum bestz — 0.0025); the third and fourth give the absolute SDSS 
2-band magnitude range; and the remaining columns contain the num- 
ber of X-ray-dctccted and X-ray-observed, optical galaxies. 

between a galaxy's spect ral energy distribu t ion a nd a 
set of template spectra. iBlanton fc Roweii (j2007[) de- 
rive five global templates from 450 instantaneous bursts 
of star formation from Bruzual fc Chario t (2003 ) spec- 
tral synthesis models, using the lChabrieil (|200l stellar 
initial mass function and the Padova 19 94 isochrones, 
and 35 templates from MAPPINGS-III (|Kewlev et alj 
[2OOI models of emission from ionized gas. These 
global templates are characterized using "nonnegative 
matrix factorization," which determines a set of non- 
negative templates that can be combined linearly to 
explain a data set. The five global templates corre- 
spond roughly to a very old galaxy template spectrum, 
a very you ng template, and several intermediate tem- 
plates (see 'Blant on fc RoweisI 120071 Figure 4). The 
kcorrect template fits can be interpreted physically 
and we utilize them to determine mass estimates for 
each of our galaxies. We sum the fit coefficients (units: 
1 Mq/{Dl/10 pc)^, where is the luminosity distance) 
and multiply by {Dl/10 pc)^ to find the total current 
stellar mass for each galaxy. In the left panel of Fig. ^ 
we plot log Mi,, in units of solar mass, as a function of 
bestz; on the right we show rest-frame (u— r)o vs. log M^,. 
To compare with studi es that use mass to define th eir 
galaxy samples (e.g., [Silverman et al.l 120081 r2009allbl ). 
we create complete mass samples, denoted samples 1* - 
5*, analogous to those in Fig. |S1 9 < log Af* < 9.8, 9.8 < 
log < 10.4, 10.4 < log < 11, 11 < log Mi, < 11.6, 
and log Mi, > 11.6 (the redshift intervals are unchanged). 
These mass intervals are selected based on a crude cor- 
respondence between Mi and log M^, i.e. one magnitude 
in Mi is roughly equivalent to 0.4 dex in log (see also 
Fig. [H. 

5. THE FRACTION 

5.1. The X-ray-Active Fraction 

We calculate the fraction of X-ray-active galaxies for 
each of the five samples complete in bestz and Mi de- 
scribed in In addition to the optical criteria, we 
apply X-ray selection criteria so that our samples are 
complete in both the optical and in X-rays. For each 
fraction, we define an X-ray threshold, L^, which gives 
the brightest X-ray luminosity to which all objects in 
the sample are, or would have been, detected by Chan- 
dra. Any detected source in the sample must have an 
X-ray luminosity at or above the X-ray luminosity limit 



(at the 90% confidence level). 

The details of the X-ray criteria are: (1) a galaxy ap- 
pears in the denominator only if its broad-band X-ray 
luminosity upper limit Lx,iim (see M.ll for a definition) 
is less than the chosen threshold, L^, regardless of X- 
ray-detection status. In other words, each optical galaxy 
counted in the denominator would have been detected 
at L'x since L'^ > Lx,iim, where Lx,iini = 47rD| fx,iim- 
(2) The numerator is the X-ray-detected subset of the de- 
nominator, where the broad-band X-ray luminosity, Lx, 
for each detection is greater than or equal to the same 
limiting L'x used to define the denominator. (We dis- 
cuss the objects that appear in the volume- limited optical 
samples, but do not pass the X-ray cut in the Appendix.) 

We can express these criteria mathematically with 
two functions, (1) G{Lx,iim, Mi, bestz) to describe 
galaxies with X-ray luminosity limits, and (2) 
X{Lx, Lx,iim, Mi, bestz) to describe the X-ray detec- 
tions. We write the denominator and the numerator of 
the fraction as: 

Nurr.= Y^ X] ^G{Lx,Um,M„bestz) , (3) 
Ndet=Yl J2 ^^X{Lx,Lx,Hm,M„bestz) , 

M/„.„t L'^ 

(4) 

where z^in, z^ax, M faint, and Mbright, define our 
volume-limited optical samples, and L'-^ establishes X- 
ray completeness. All five volume- limited samples have 
a minimum redshift Zmin = 0.0025. 

Figure [10] shows the X-ray-active fraction as a function 
of L^, where the luminosity threshold is set to the mea- 
sured Lx for each of the X-ray-detected sources; panels 
1-5 show the fraction for the five volume- limited opti- 
cal samples. Errors on the fraction due to small number 
statistics, estimated via the /3 distribution (gray filled 
circles; see §5.21 below), are largest at the lowest lumi- 
nosity limits, where both the number of X-ray detections 
and galaxies is small (see also Fig. [TH]in the Appendix). 
The sixth panel shows the active fraction over a narrower 
range of L'-^ from all five samples on the same vertical 
scale (sample 1: blue solid line; sample 2: green dotted 
line, sample 3: purple dashed line; sample 4: magenta 
dot-dashed line; sample 5: red triple dot-dashed line). 
The dashed gray line denotes an X-ray luminosity limit 
of L'j^ = lO''^ ergs~^, above which we consider our frac- 
tion a proper AGN fraction (see and below which 
the fraction may be impacted by a mix of stellar emis- 
sion (principally from X-ray binaries) and hot ISM in the 
host galaxies. 

The redshift distributions of the five optically-complete 
samples differ substantially from one another, see Fig. 
[TTl Sample 1 (at the lowest rcdshifts) shows two peaks, 
one corresponding to a large number of spectroscopic red- 
shifts at z = 0.02 and the other dX z = 0.075; sample 2 
peaks at z = 0.23; sample 3 at z = 0.35; sample 4 at 
z = 0.46; and sample 5 at about z = 0.6. As a result, 
the increase in X-ray-active fraction between samples 1 
and 5 may be due either to a trend with Mi, with z, or 
with both. We discuss these dependencies in i)5.3l and ^ 
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Fig. 8. — Absolute i-band magnitude vs. bestz for the 105,550 
galaxies in the X-ray-observed sample (light gray points) and the 
1,593 X-ray-detected sources (large dark gray stars). Colored lines 
and large symbols / nu mbe rs indicate the five independent, volume- 
limited samples (see i|4.2l and Table [2]l for which we determine the 
X-ray-active and AGN fractions. Representative 1 cr error bars are 
as described in Fig. |6] 

(A color version of this figure is available in the online journal.) 



5.2. Error Estimation via the j3 Distribution 

We face small number statistics throughout this study, 
particularly when we bin our modest number of X-ray- 
detections further by mass, redshift, color, etc. Statistics 
used to estimate errors on numbers of ev ents, e.g., the 
Poisson distribution or the ' Gehrelsl (jlQSGT ) statistic, are 
not suitable for estimating errors on fractions. When 
it comes to estimating uncertainties on a fraction, one 
should use the binomial distribution, or its conjugate 
prior, the (3 distribution (jEvans et al.l l200"ot ) . Given a 
fraction, F , and a number of trials, iV, the binomial dis- 
tribution describes the probability of different outcomes. 
Given the outcome of such a trial, the /3 distribution de- 
scribes the probability distribution for the fraction itself. 
Hence, we determine both the X-ray-active (or AGN) 
fraction and the uncertainties on the fraction using the 
regularized incomplete (3 distribution (see Equation [S] 
and Fig. [pi)- 

The /3 distribution is continuous, yields a fraction that 
is well behaved in the extremes where either the numer- 
ator or the denominator (or both) are zero, and properly 
accounts for a numerator and denominator that are not 
drawn from independent samples, i.e. where one is a sub- 
set of the other, as is the case here. We assign a prob- 
ability distribution to the fraction, (defined on the 
interval [0, 1]), given the measured values a = Ndet + 1 



and /3 = Nu. 



N, 



det ' 



1. For each input Ndet and Nu 



we calculate a cumulative distribution function (CDF) 
via the /3 distribution and choose the mean of the CDF 



TABLE 3 
X-RAY Active Fraction 



log L' 


^det 




^ det, drop 


^ lim,drop 


min 


mean 


max 


Sample 1 


-18 > M, > 


-20, bestz < 0.125 








41.0 


7 


1587 


7 


3441 


0.33 


0.50 


0.68 


41.5 


8 


2974 


3 


2054 


0.20 


0.30 


0.40 


42.0 


6 


4336 





692 


0.10 


0.16 


0.22 


42.5 


1 


4930 





98 


0.01 


0.04 


0.07 


43.0 





5026 





2 


0.00 


0.02 


0.04 


43.5 





5028 








0.00 


0.02 


0.04 


44.0 





5028 








0.00 


0.02 


0.04 


Sample 2 


-20 > > 


-21, bestz < 0.275 








41.0 


2 


420 


77 


10529 


0.33 


0.71 


1.10 


41.5 


11 


1629 


57 


9320 


0.53 


0.74 


0.94 


42.0 


21 


4415 


30 


6534 


0.39 


0.50 


0.60 


42.5 


24 


7930 


7 


3019 


0.25 


0.32 


0.38 


43.0 


7 


10240 





709 


0.05 


0.08 


0.11 


43.5 





10903 





46 


0.00 


0.01 


0.02 


44.0 
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Sample 3 


-21 > M, > 


-22, bestz < 0.400 
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1.96 


2.36 
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Sample 4 


-22 > Ah > 


-23, bestz < 0.550 
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20.45 
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2.61 


3.53 


4.44 
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2.47 
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3.24 
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1.49 


1.65 


1.81 
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44 
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0.59 


0.66 


0.73 


43.5 


44 


19154 
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2941 


0.20 


0.23 


0.27 


44.0 
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21732 
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0.01 


0.02 


0.03 


Sample 5 


M, < 


-23, bestz < 0.700 
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18 
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49.10 


60.00 


70.87 


41.5 


13 


84 


291 


13544 


12.34 


16.28 


20.23 


42.0 


15 


419 


260 


13209 


2.88 


3.80 


4.72 


42.5 


42 


1842 


199 
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1.98 


2.33 


2.68 


43.0 


96 


5499 


88 


8129 


1.59 


1.76 


1.94 


43.5 


94 


10085 


15 


3543 


0.85 


0.94 


1.04 


44.0 


51 


12875 


2 


753 


0.35 


0.40 


0.46 



Note. — The first column gives the broad-band (0.5-8.0 kcV) X-ray luminos- 
ity threshold. The second and third columns contain the number of detections 
and limits that were used to calculate the mean X-ray-active fraction, Fmean', 
columns four and five contain the number of objects that did not pass the X-ray 
luminosity cut (details in the Appendix). Fmin and Fj^ax are the 1 — cr upper 
and lower limits on the X-ray-active fraction (sec jJOJ. The bold rows indicate 
the AGN fraction, i.e. the fraction for log L ^ — 42. 

as the best value for the X-ray-active fraction, such that 

a Ndet + 1 



F„ 



a + /3 



(5) 



The asymmetric 1 a error^ are calculated based on the 
median value of F (evaluated where = 0.5), and de- 
fined such that Fmax corresponds to Pp{0.5) + 1 a and 
Fmin corresponds to Pg (0.5) — 1 cr (see Fig. [T2]for an 
example). Throughout the remainder of this manuscript 
(tables excepted), we specify the asymmetric 1 a errors 
only in cases where the upper and lower errors differ by 
more than 10%; in all other cases we quote the larger 
of the two. (That is, if we have F±yj "^iH quote 
F ± max{X,Y), unless X and Y differ by more than 
10%.) The cumulative X-ray-active fractions and their (3 
distribution errors, for a set of seven useful X-ray lumi- 
nosity thresholds, appear in Table [3l 
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Fig. 9. — (Ijeft) Log of the host galaxy mass, log M*, vs. bes tz (g ray symbols as in Fig. |8]l. A set of samples, complete in galaxy mass, are 
indicated with colored lines and large symbols/numbers (see i|4.2l l. These arc analogous to the volume-limited samples in Fig. \E\ (Right) 
Rest-frame [u — r)o color vs. log Af* (gray and black symbols as in Fig. |6ll. Black lines mark log Af* = 10.6 (solid) and 10.4 (dashed) and 
are used to define mass limits for a comparison between our field fractions and those of other X-ray surveys (see Q. Representative 1 cr 
error bars in both panels are as described in Fig. [6] 
(A color version of this figure is available in the online journal.) 



5.3. The AGN Fraction 

In this section and the discussion that follows, we dif- 
ferentiate between the X-ray-active fraction, which can 
be defined for an arbitrary threshold, and the AGN 
fraction, Fagn- For our purposes, the AGN fraction will 
refer only to the cases where the limiting X-ray threshold 
is L 'y = 10"^^ ergsec" ^, a typical cutoff for AGN activ- 
ity (Zezas et al.lll998| ). For each of our volume-limited 
samples, we measure the following AGN fractions (see 
also Table H bold entries): Fagn i = 0.16 ± 0.06%, 
Fagn.2 = 0.50 ± 0.11%, Fagn.3 = 1.27 ± 0.18%, 
Fagn A = 2.85 ± 0.39%, and Fagn^s = 3.80 ± 0.92%. 

To test for evolution in the AGN fraction with bestz, 
and/or a change with Mi or galaxy mass, we further sub- 
divide our five volume-limited samples, re-calculate the 
AGN fraction in these finer bins, and perform a Spear- 
man rank correlation analysis (including errors) to as- 
sess the significance of any trends. We first sort our data 
into four absolute magnitude bins, —20 > Mi > —21, 
-21 > M,> -22, -22 > M, > -23, -23 > M, > -24, 
and simultaneously fit all twelve of the data points shown 
in the left-most panel of Fig. [T31 We find a marginally 
significant trend at the ^76% confidence level between 
the AGN fraction and bestz. In four redshift bins, 
z < 0.125, 0.125 < z < 0.275, 0.275 < z < 0.4, and 
0.4 < z < 0.55, we find that the AGN fraction correlates 
strongly with Mi, at the 91% confidence level (96% if the 
0.4 < z < 0.55 redshift bin is ignored; Fig. [T21 middle 
panel). For the same redshift bins, the correlation with 
galaxy mass is similar, but slightly stronger, ~93% and 
96% with and without the higliest-redshift bin (Fig. [T3l 
right panel). 

These trends reveal a field AGN fraction that increases 



toward brighter absolute magnitudes and higher masses, 
as might be expected from cosmological models wherein 
the merger of massive galaxies fuels both accretion onto 
a central SMBH and star formation. The weaker evolu- 
tion of the AGN fraction with redshift indicates that the 
mechanism for black hole fueling may be changing, e.g., 
shifting from high luminosity merger-driven activity to 
fueling via tidal interactions, bars, d ynamical frict ion, or 
other lower luminosity drivers (see iMartinil 12004 for a 
review). In such a scenario, the lower AGN fraction may 
reflect the evolution of the AGN X-ray luminosity func- 
tion, which is characterized by a decrease in accretion lu- 
minosity toward the present day. In the left-most panel 
of Fig. [T3]we show two common fits to the AGN X-ray lu- 
minosity function (lUeda et al.ll200l iBarger et all 120051: 
iHasinger et"alll2005r normalized to the mean AGN frac- 
tion in the four Mi bins at z ~ 0.21: pure luminosity 
evolution cx (1 + z)'^ (gray dotted line), and luminosity- 
dependent density evolution cx (1 -I- z)* (dashed gray 
line). We explore the possibility of an evolving AGN 
fraction further in ^J6l 

5.4. Redshift Monte Carlo Simulations 

As discussed in %.2l small number statistics for X-ray- 
detected sources in our volume-limited bins place a basic 
limit on our ability to accurately determine the X-ray- 
active fraction. However, our use of photometric red- 
shifts for the majority of the galaxies (^ 97%) in our X- 
ray-observed sample introduces another potentially large 
source of error. Photometric redshift errors are of or- 
der 3%, but cai i be as high as 10% at fainter appar- 
ent magnitudes (jCsabai et al.ll2003l : lOvaizu et aLl HoOS). 
Meanwhile, photometric redshifts for quasars can experi- 
ence catastrophic failures due to aliasing ((Budavari et 
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Fig. 10. — The first five panels show the X-ray-active fraction as a function of limiting broad-band X-ray luminosity for each of the five 
volume-limited samples described in Table [2] and shown in Fig. [S] The solid lines connect the fraction as determined for the limiting L'^ 
at each X-ray detection in the sample. The 1 a errors (determined via the (3 distribution; see i|5.2l l appear at each point where an active 
fraction was calculated. The final panel contains the best value for the fraction for each sample, plotted on the same scale for ease of 
comparison; sample 1 (solid blue line), sample 2 (dotted green line), sample 3 (dashed purple line), sample 4 (dot-dashed magenta line), 
and sample 5 (triple dot-dashed red line). The dashed gray line in this last panel denotes the limiting X-ray luminosity, log Lx = 42, that 
we associate with AGN; at this limit the AGN fractions are Fagn,i = 0.16 ± 0.06%, Fagn,2 = 0.50 ± 0.11%, Fagn,3 = 1-27 ± 0.18%, 
Fagna = 2.85 ± 0.39%, and Fagn,5 = 3.80 ± 0.92% for samples 1-5, respectively (see also Table^. 
(A color version of this figure is available in the online journal.) 
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Fig. 11. — Distribution of bestz for each of the five samples de- 
scribed in Table [2] and in Fig. |8](line styles are as in the bottom, 
right panel of Fig. IIOI I. 

(A color version of this figure is available in the online journal.) 




Na., = 50, 
N„„ = 20000. 
F„;„ = 0.2195 
= 0.2533 
= 0.2905 



0.20 0.25 

X-ray-actlve Fraction [%] 

Fig. 12. — The cumulative f} distribution function (black solid 
line) for a toy example where the number of detections, Ndet, 
equals 50 and the number of X-ray observed galaxies, the denomi- 
nator Nii^, is 20,000. For these inputs, the mean fraction is 0.25%. 
The mean value of the CDF, i^meam is marked with a large, filled 
circle, as are the slightly asymmetric 1 cr errors (also marked with 
dashed lines), corresponding to F^in and Fmax- 



120011: iRichards et al.|[200l . Hence, photometric redshift 
errors dominate the uncertainty on our volume-hmited 
samples. To estimate the impact of the photo-z uncer- 
tainties on our active fraction, we perform two Monte 
Carlo (MC) simulations: a "full" MC and a "quick" MC. 
Pull MC: In the full Monte Carlo simulation, the 1 cr 
error for each object not associated with a QSO pho- 
tometric redshift is multiplied by a normally-distributed 
random number selected from a distribution with a mean 
of zero and a standard deviation of one. This quantity is 
added to the best redshift, yielding a new hestz. 
Since single- valued quasar photometric redshifts, de- 



rived using broad band filters, suffer from aliasing in 
certain redshift regimes they are better described by 
their full redshift probabilit y distribution fu nction (see, 
e.g.. IMvers et al.| [2009: Ric hards et all 120091 and refer- 
ences therein). Thus, for objects where we have as- 
signed a QSO photometric redshift we use the full red- 
sh ift PDF, obtained using the nearest neighbor approach 
of iBall et al.l ()2008[ ). to calculate the new hestz. We also 
run the MC tests using the QSO 1 a errors described 
in 3321 instead of the full PDF. The results of those 
tests show excellent agreement with the MCs that use 
the QSO PDF. Thus, we do not attempt to acquire full 
redshift PDFs for all of the objects in our sample. 

We use the new hestz values for non-QSO and 
QSO objects to determine new luminosity distances, k- 
corrections, and absolute magnitudes via kcorrect. We 
then pull five volume-limited samples analogous to those 
used for our primary analysis using the new bestz and its 
associated Mi. We adopt the same seven X-ray luminos- 
ity thresholds and calculate new X-ray-active fractions 
for each. 

The above constitutes a single iteration of the full MC. 
The process is computationally expensive; determining 
the spectral template match and k-correction for each 
of ~ 100,000 objects is the largest computational task. 
Thus, we perform only 100 iterations of the full MC, 
which result in 100 possible X-ray-active fractions for 
each volume-limited sample, at each of seven X-ray lu- 
minosity limits. We take the median and standard de- 
viation of these fractions and use them to estimate the 
uncertainty contributed by the photo-z errors. We also 
use the full MC to estimate the representative error bars 
shown in Figs. 6-9. 

Quick MC: In our quick Monte Carlo simulation, new 
bestz values are derived as in the full MC and new lu- 
minosity distances are then calculated. We assume that 
a small change in the redshift will result in only a small 
change in the k-correction and so leave the k-correction 
unaltered. Hence, the new absolute magnitude is simply 
the original Mi weighted by the ratio of the luminosity 
distances: M^^^ew = - 5 log [DL^new/DL). We per- 
form this quick MC initially for 100 iterations and make a 
direct comparison to our full MC. We then perform 1000 
iterations to achieve more robust statistics. Here again, 
we derive 100 (or 1000) X-ray-active fractions and com- 
pute a median fraction and standard deviation for each 
volume-limited sample at each of seven X-ray luminosity 
thresholds. 

The errors that result from the redshift MCs are, in 
general, smaller than the /? distribution errors, though in 
some cases they are of comparable magnitude. For exam- 
ple, if we adopt an X-ray threshold of L'^ — 10'*^ ergs"^ 
for sample 3, the 100-iteration full MC gives a median 
fraction of 1.25 ± 0.12%, the 100-iteration quick MC 
yields a median of 1.25 ± 0.11%, and the 1000-itcration 
quick MC gives a median fraction of 1.24±0.11%. These 
compare to a /? distribution result of 1.27±0.18%, where 
the 13 distribution errors arise from small numbers of X- 
ray detected galaxies, instead of photo-z uncertainties. 
The fi distribution and MC errors could be combined in 
quadrature; or as a conservative estimate, the errors in 
Table [3] might simply be increased by a factor of V2 to 
take photometric redshift errors into account. 
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Fig. 13. — AGN fraction, i.e. L'^ = 10*^ ergsec in bins of redshift (left), absolute i-band magnitude (center), and mass (right). Upper 
limits on the fractions are plotted as solid, downward-facing black triangles (see ^5.51 for details). Two common fits for the evolution of 
the X-ray luminosity function of X-ray-detected AGN, normalized to the mean AGN fraction in the four Mi bins at 2: ~ 0.21, are shown 
in the left-hand panel; pure luminosity evolution oc (1 -|- z)"^ (gray dotted line), and luminosity-dependent density evolution <x (1 -|- 2)^ 
(dashed gray line). The Spearman rank correlation test indicates that the active fraction correlates most strongly with galaxy mass (93% 
confidence; 96% if the 0.4 < 2 < 0.55 redshift bin is ignored). There is also a correlation with Mi (91% and 96% confidence, respectively). 
A weaker correlation with bestz is detected at only the 76% confidence level. (A color version of this figure is available in the online journal.) 



5.5. Selection Effects and Sample Bias 

As discussed in ^3.21 we have included a large num- 
ber, (~27 thousand) of red, [u — g) > 0.6, point sources 
in our galaxy sample to insure that we do not bias our 
AGN fraction at faint apparent magnitudes (and at high- 
redshift) where morphological criteria are less reliable, 
and to include all possible QSOs. As a result, the de- 
nominator of our fractions is likely to be contaminated 
by stars and the AGN fractions reported here must be 
taken as a lower limit to the "true" AGN fraction. 

To quantify the magnitude of the stellar contam- 
ination, we return t o SDSS star-galaxy separation 
(jStoughton et al.ll2002f ). We assume that all of the point 
sources with red colors (u — g > 0.6) are actually stars 
and remove these 27,958 objects along with their X- 
ray counterparts (120 detections) from our galaxy sam- 
ple. We re-run our analysis and find mean AGN frac- 
tions at L'x = 10^*2 ergsec-i of: 0.16 ± 0.08% (sample 
1), 0.61 ± 0.14% (sample 2), 1.53 ± 0.33% (sample 3), 
3.45 ± 0.46% (sample 4), and 5.71 ± 1.37% (sample 5). 
Since some of these excluded objects are certainly galax- 
ies, and the number hosting X-ray-emitting AGN is likely 
to be smaller than the passive population, we consider 
these fractions to be upper limits, i.e. we will have elim- 
inated more legitimate galactic objects from the denom- 
inator than the numerator. We plot these "maximum" 
AGN fractions in Fig. [13] (filled, black, downward-facing 
triangles) and note that, with few exceptions, they agree 
with our original AGN fractions within the errors. 

We also consider the possibility that the combination of 
ChaMP and SDSS spectroscopy (which includes known 
X-ray sources in its targeting algorithms) introduces a 
bias in our redshifts, leading to greater spectroscopic 
completeness in the X-ray-detected sample than in the 
X-ray-observed sample. This is, at worst, a second- 



order effect impacting only the precision of the redshifts 
since we include galaxies in our sample whether or not 
they have spectroscopy. One exception is along the stel- 
lar locus. In regions of the color-color diagram where 
we have eliminated objects based on their photomet- 
ric redshifts and colors (see §5.5), we do not eliminate 
spectroscopically-confirmed galaxies. Thus, we may have 
an over-representation of objects in these particular color 
regimes due to existing ChaMP and SDSS spectroscopy. 

Pollution of the host galaxy light by a central AGN 
or active star formation is another source of bias in 
our galaxy sample. Though we are interested in locat- 
ing active black holes in nearby galaxies, these same 
sources may contaminate the spectral energy distribu- 
tion (SED) of the host galaxy and thus impact the pho- 
tometric redshift algorithms and kcorrect spectral tem- 
plate matching. There are several possible approaches 
to this problem: (1) in principle it is possible to correct 
the host galaxy SED by extrapolating an AGN power- 
law spectrum, e.g., normalized to the measured X-ray 
luminosity, into the optical and then subtracting this 
component. This method is not likely to be efficient 
for our sample because only a small number (^ 1.5%) 
of our galaxies have X-ray detections. Also, a non- 
negligible fractio n of X-ray AGN have the optical SEDs of 
normal galaxies (Georgantopoulos & Georgakakis"2005t 
Hornschcmcicr et al. 2005; Kim ct al. 2006; Civano ct al] 
(20071 : iCocchia et al.ll2007 ). so subtraction of a canonical 
normalized SED could introduce other errors. (2) Where 
spectra are available, the AGN contribution can be es- 
timated and removed from the host galaxy SED. With 
limited spectroscopic coverage (~ 2.5% of our galaxies 
have spectra) , this too is an imperfect approach for these 
data. (3) Lastly, fitting of high-quality imaging data can 
successfully deconvolve the AGN point source from the 
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host galaxy light. However, since the large majority of 
our objects are faint such an approach is not feasible. 
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Fig. 14. — Absolute i-band magnitude vs. log Af* for the lim- 
its (small points) and detections (stars) in the Full galaxy sample. 
Th e sym bols are colored according to the {u — r)o criteria described 
in i|6.1l The strong correlation between magnitude and mass is ex- 
pected as more massive galaxies tend to he more luminous. In 
addition, both are calculated here from the fit performed by the 
kcorrect tool jBlanton &: Rowei£ii2007i ). There is a population of 
objects that define a "second" locus or whose Mi or mass determi- 
nations are flawed. We differentiate these objects by eye (dashed 
line) and consider the possibility that thes e con tain a central AGN 
that dominates the host galaxy SED (see il5.5ll . 

We instead make a qualitative assessment of how an 
AGN contribution to our host galaxy SED would bias our 
fractions. In Fig. [T31 we plot the i-band absolute magni- 
tude versus the stellar mass of our galaxy sample (small 
points) and the X-ray-detected subset (stars), color- 
coded according to the following {u — r)o color sequences: 
blue: {u — r)o < 1.8, green: 1.8 < (u — r)o < 2.6, and 
red: (u — r)o > 2.6. Since more massive galaxies tend to 
be more luminous, it is not remarkable that these quanti- 
ties show a strong correlation. In our analysis, the galaxy 
magnitudes and masses are de rived using the sarne fit to 
the models implemented in the lBlanton fc RoweisI (|2007[ ) 
kcorrect tool, further insuring a correlation between the 
two. We note a distinct second population of blue ob- 
jects with apparently higher masses and brighter lumi- 
nosities; we demarcate this "second sequence" by eye (as 
indicated by the black dashed line in Fig. [T31). Using 
this simple cut, we identify 3,023 X-ray-observed and 
156 X-ray-detected galaxies where AGN contamination 
may be an issue. These objects are indicated in Figs. 
[6] and [14] (black points [observed] and black open cir- 
cles [detected]). Where we have only upper limits to the 
X-ray luminosity, the galaxies' very blue colors are our 
primary evidence that they host an AGN. To be conser- 
vative, here we assume that all of these galaxies do host 
an AGN whose X-ray emission is either weak or obscured. 



but whose optical continuum substantially contaminates 
the colors. 

We might explain this second sequence by speculating 
that the k-correction and derived Mi and stellar mass 
for these objects is incorrect. For example, we might 
estimate, based on Fig. [SJ that the derived Mi for these 
objects is too faint by about one magnitude. Hence, for 
an object at Mi ^ —23, we would correct up by one 
magnitude to Mi ^ —24. An object on the primary 
sequence at this magnitude corresponds to a stellar mass 
of log Mstar ^ 11.8, very reasonable for an AGN host. 

Mischaracterization of the k-correction is also likely to 
bias our AGN fraction as a function of color, yielding an 
artificially high fraction in the blue sequence. We have, 
however, tested the impact of these "second sequence" 
objects on our global AGN fraction as a function of Mi 
and hestz and found that the fraction is robust to their 
inclusion (or exclusion). 

Among these potentially contaminated objects, we as- 
sociate 121 (4%) with QSO photometric redshifts — this 
compares to ^ 0.3% of objects with QSO photometric 
redshifts in the Full sample. Another 83 have spectra, of 
which 66 have secure spectroscopic types: there are 50 
broad line AGN (BLAGN) identified via ChaMP spectro- 
scopic follow-up, a single ChaMP narrow emission line 
galaxy, and 15 SDSS MPA/JHU narrow emission line 
galaxies. Thus, a significant fraction of the objects along 
this second locus are indeed associated with active, star- 
forming galaxies and/or likely AGN. However, since these 
objects lie preferentially outside of the limits of our five 
volume-limited samples — only one X-ray-detection and 
fewer than 300 galaxies across all five samples appear 
in our AGN fraction calculations — we do not recalcu- 
late Fagn without these objects. We do, however, cau- 
tion that AGN activity may bias our derived host galaxy 
properties in a small number of cases and in the highest 
redshift, highest absolute magnitude sample, in particu- 
lar. 

White dwarfs are another possible contaminant in our 
galaxy sample. Their blue colors and X-ray brightness 
can cause them to masquerade as AGN, particularly if 
they have been incorrectly assigned a cosmological pho- 
tometric redshift. However, since we removed objects 
with significant proper motions, we anticipate that these 
will not contribute more than an object or two to our 
numerator and thus will not seriously bias our fractions. 

6. DISCUSSION 

In this section we compare our results with previ- 
ous AGN fraction measurements in the field and in 
groups/clusters. We also expand our analysis of the AGN 
fraction as a function of host galaxy properties. To ex- 
tract Vega UBVRI absolute magnitudes, k-corrections, 
and masses, we re-run the kcorrect tool on our Fu ll 
galaxy sample (|Besselll ITQQOI : iBlanton fc Roweid l2007[ ) : 
these photometric bands are useful for comparisons with 
other studies in the literature. 

6.1. Comparison to Previous Field Fractions 

iPagd (j2001h combines the present-day optical luminos- 
ity function of galaxies with the X-ray luminosity func- 
tion of Seyfert Is to calculate an expected AGN fraction. 
To reproduce the local type 1 AGN X-ray luminosity 
function for a plausible mix of AGN hosts, including both 
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early- and late-type galaxies, he finds a Seyfert 1 fraction 
of about 1% at 2; ~ 0. This compares favorably to our di- 
rect measurement using intermediate-luminosity galaxy 
sample 3 (—21 > Mi > —22) at low-redshift, for which 
we find an AGN fraction of 1.27 ± 0.18%. 

From a compilation of deep X-ray surveys (including 
the Chandra Deep Fields, the Chandra Large- Area Syn- 
optic X-Ray Survey of Lockman Hole [CLASXS] , and the 
All-wavelength Extended Groth strip International Sur- 
vey [AEGIS]; covering - 2 deg^l. lShi etall ^0§) derive 
an X-ray AGN fraction as a function of host stellar mass. 
We present the AGN fraction versus host mass in Fig. 
[15] for our ChaMP/SDSS sample (color ed symbols indi - 
cate our five redshift intervals) and from lShi et al.l ()2008f ) 
(black symbols; crosses for 0.1 < z < 0.4 and diamonds 
for 0.4 < z < 0.7). The fractions from these studies agree 
to within 1 cr in all but one case, the highest mass bin of 
their low-redshift sample, where the agreement is good to 
within 1.5 a. We also include a comparison to the power- 
ful AGN fraction from an optic al, spectroscopic study o f 
SPSS galaxies ( solid black line; IKauffmann et al]|2003[ ). 
I Shi et afl ([2008) re-normalize this relation by a factor of 
8 to fit their data and account for differences between 
the X-ray and optical selection functions, e.g., emission- 
line-selected AGN samples are more complete (less sensi- 
tive to obscuration) for low- mass galaxies at low-r edshift 
than X-ray-selcctcd AGN ( Heckman et al.ll2005D . The 
[Kauffmann ct al. (2003) SDSS powerful AGN fraction is 
measured for 2: < 0.3 and shows an increase with stellar 
mass to 10^ and then a sharp de crease toward galax- 
ies of higher masses. iShi et al] (|2008f ) interpret their low- 
redshift sample as demonstrating a turn over (or flatten- 
ing) in the AGN fraction near IO^^A'Iq, though no such 
inflection is detected for their high-redshift sample. Nei- 
ther our low- or high-redshift samples show a turn over 
in the fraction, but instead increase from the lowest-mass 
to the most massive galaxies in our study. 

We also compare o ur field AGN fraction t o the zCOS- 
MOS survey results ([Silverman et al.ll2009 aV which uti- 
hze XMM-Newton observations of the COSMOS fields. 
For a mass limit of log M^, > 10.4 (black dashed line in 
the right panel of Fig. [9]) we m atch the redshift and X- 
ray luminosity intervals used bv [Silverman et ahl (|2009afl 
in their Fig. 11: 0.1 < z < 0.58 and 0.58 < z < 1.05, 
with Lx > 10^^'^ ergs~^ (dashed black line in the right 
panel of Fig. [9]). In the low redshift bin we have 186 X- 
ray detections and 19,057 X-ray-observed galaxies; these 
yield an AGN fraction of 0.98 ± 0.07%. The zCOSMOS 
field fraction in this bin is 1.9 ± 0.6%. Hence, our field 
AGN fraction faUs below the zCOSMOS field fraction by 
1.5 a. The difference between the ChaMP/SDSS and 
zCOSMOS low-redshift results may reflect evolution in 
the AGN fraction if the zCOSMOS sample is weighted to- 
ward higher-redshift objects (in the low-redshift bin the 
ChaMP/SDSS sample has an average redshift of (z) — 
0.28). In the high redshift bin, we flnd 14 detections and 
296 limits, which yield a fraction of 5.03 ± 1.25%, vs. the 
zCOSMOS value of 3.7±g:i%. However, we do not con- 
sider this second a robust comparison because our sample 
is not complete beyond redshifts of 0.4 for this mass cut 
(see Fig. [HI left panel). 

Evidence for evolution i n the AGN fraction h as been 
seen in other recent work. iLehmer et al.l ([2007[ ) use the 
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Fig. 15. — Log of the AGN fraction (log L'^ = 42, solid black line 
in the right panel of Fig. [6]l versus stellar mass. Colored symbols 
indicate our AGN fractions for five redshift intervals {z < 0.125, 
blue filled circles; 0.125 < 2 < 0.275, green open circles; 0.275 < 
z < 0.4, purple filled squares; 0.4 < 2 < 0.55, magenta open 
squares; 0.55 < 2 < 0.7, red filled triangles) from the mass-selected 
sam ples of Fi g. [9] (left panel). The high- and low-redshift fractions 
from lShi et~al. (20Q8 ) are plotted as black symbols (0.1 < 2 < 0.4 
as crosses, 0.4 < 2 < 0.7 as open diamonds). The solid line shows 
the trend of the SDSS powerful AG N fraction with the host stellar 
mass from Kauff mann et al.l ( |2003l ) ; the dashed line shows the same 
trend wi th the normalization decreased by a factor of 8 (following 
IShi et al.ll2008i ). 

(A color version of this figure is available in the online journal.) 

Chandra Deep Fields to study the AGN fraction in early- 
type galaxies (to z 0.7) and find evolution consis- 
tent with the (1 -I- z)^ pure luminosity evolution (PLE) 
model frequently fit to the luminosity functi on of X-ray- 
selectcd A GNs (hereafter the AGN XLF: lUeda et"all 
[2003: Barger et al.|[20Ql iHasinger et al.ll2005[ ). Similar 
behavior is also noted in an earlier analysis of the stacked 
X-ray properties of early-type galaxies in a 1.4 deq'^ 
field in the NOAO Deep Wide-Field Survey ([Brand et al.l 
'2005") . (Note, however, in a study of late-type galaxies, 
Lehm er et al.l ((2008) find no evidence for evolution of the 
AGN fraction over the redshift range 0.1 < z < 0.8.) In 
the left panel of Fig. [T3| we show curves for PLE (cx 
(1 + z)^; gray dotted line) and for luminosity-dependent 
density evolution (LDDE), another popular model for 
the evolution of the AGN XLF (oc (1 4- z)^; gray dashed 
line). These curves are intended to be illustrative only 
and are thus arbitrarily normalized to the mean AGN 
fraction in the four Mi bins with 0.125 < z < 0.275. For 
the —22.0 > Mi > —23.0 magnitude bin, we compare the 
AGN fraction in the lowest and highest redshift intervals 
and find, Fagn{z = 0.47) = 2.2±t?-5 xFAGNiz = 0.09); 
i.e. though poorly constrained, the change in this AGN 
fraction is consistent with either (1 + z)^ or (1 + z)^ evo- 
lution. 

In the current work, we predicate our comparison to 
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TABLE 4 
Comparison to Cluster Fractions 



Mji < min max N^j^^ N;j„ min mean max ^clust C^) 

-20 40.9 43.6 26 448 4.89 6.00 7.11 6.0 

-20 41.0 ... 28 642 3.69 4.50 5.31 4.9 

-20 42.0 ... 113 9588 1.08 1.19 1.30 1.0 

-21.3 41.0 ... 25 208 10.12 12.38 14.64 9.8 

-21.3 41.0 ... 4 46 6.14 10.42 14.71 5.5 



Note. — " The field fractions calculated from the ChaMP and SDSS survey 
data and described in the present work for 0.0 5 < 2 < 0.31. The first four 
cluster fractions come fromlM^rt^^et (2007) and also span 0.05 < z <. 0.31; 
the final comparison is to [^m^^^n ~j2007) clusters with 0.01 < 2: < 0.05. 
Absolute R magnitudes arc calculated via the SDSS transformations described 
in ^4.11 The Sun ct al. (2007) result in the table is a simple fraction determined 
from their values of 9 AGN and 163 galaxies. The mean AGN fraction for these 
inputs, calculated via the beta distribution, is 6.1%. 



the AGN XLF on the assumption that the underlying 
galaxy population does not evolve substantially between 
redshifts of z = 0.5 and the present day — indeed, we as- 
sume no evolution in our galaxy k-corrections. Thus, the 
evolution of the AGN fraction is entirely attributed to the 
evolution of the X-ray-detected AGN population. The 
assumption of zero evolution for the galaxy population 
is almost certainly too simple, particularly for late-type 
galaxies, which have been shown to fad e by a magnitude 
or more between z ~ 1 and the pres ent (jLillv et al.lll995l : 
IWolf et alJ [20031 : IFaber et all 120071) . In fact, the lack of 
evolution of t he AGN fraction in late-types reported by 
iLehmer et all f2008), may indicate that the AGN and 
late- type galaxy populations co-evolve for z < 1. 

6.2. Comparison to Cluster and Croup X-ray Fractions 

The impact of environment on the formation and evo- 
lution of accreting supermassive black holes remains an 
important open question, i.e. are AGN preferentially 
located in field galaxies, in galaxy groups, or in rich 
clusters? Early studies reveale d a paucity of optically- 
luminous AGN in clusters (Dr essier et al.l 11985) , com- 
pared to the field — 7% vs. 31% of galaxies had emission- 
line nuclei in clusters vs. the field , and 1% vs. 5% har- 
bored AGN. iGiovanelli fc HavnesI (|1985f ) proposed that 
galaxies in clusters lose cold gas from their disks by in- 
teracting with the cluster core, leaving behind a clus- 
ter galaxy population deficient in the cold gas neces- 
sary to fuel AGN. A smaller frequency of AGN fueling 
events was also proposed — in this scenario, the merger 
of two gas-rich galaxies was invoked as the mechanism 
for d riving cold inflows toward a massive central black 
hole (H arnes &: Hernauis^ll992f ). If galaxies in the clus- 
ter environment experience fewer mergers due to their 
high velocity dispersions, less AGN activity would re- 
sult. These two scenarios might also work in tandem to 
produce smaller AGN fractions in clusters. (Note that 
groups, which have lower velocity dispersions than clus- 
ters, should have merger rates closer to those in the field.) 

The last decade has shown great advances in our un- 
derstanding of the X-ray- selected populat i on of AGN 
in clu sters of galaxies iMartini et all I200a 120061 120071 
l2009f ). Enabled by the excepti onally high sen s itivity and 
spatial resolution of Chandra^ iMartini et al.l (|2002f) and 
iMartini et al.l ()2007f ) showed that the fraction of X-ray- 
selected AGN observable in clusters at 0.05 < z < 0.31 



was, in fact, higher than previously reported and simi- 
lar to the field, though the field fraction has historically 
been poorly constrained. 

To compare our robust field X-ray AGN fractions to 
those reported for clusters, we employ absolute R mag- 
nitudes, Mji^ adopt a variety of L'-^ and Mji limits, 
and measure the active fraction for 0.05 < z < 0.31 
(to m atch the cluster redshift range of IMartini et al.l 
l2007f l . These comparisons are summarized in Table HI 
Our field fractions for Mr < -20 are: (1) 6.00 ± 1.11% 
(8 X 10^° < Lx < 4 X lO'^s ergs-i); (2) 4.50 ±0.81% for 
a limiting Lx = lO^^ ergs"!; and (3) 1.19 ± 0.11% for 
a limiting Lx = lO'*^ ergs~^. These compare to cluster 
active fractions of 6.0%, 4.9%, and 1.0%, respectively. 
This excellent agreement between the field and cluster 
fractions is also borne out for a brighter magnitude cut. 
For Mfl < —21.3 and a limiting Lx = 10^^ ergs~^, we 
find a fraction of 12 .38 ± 2.26%, which agrees with the 
IMartini et cll (|2007D value of 9.8% within 1.5 a. For 
these same magnitude and X-ray luminosity limits, we 
also find a fraction of 10.42 ± 4.29% in the narr ower red- 
shift range (0.01 < z < 0.05) used bv iSun et a l. (2007), 
who find a cluster AGN fraction of 5.5%. The striking 
concordance between cluster and field fractions out to 
z ^ 0.3 implies that AGN activity at low-redshift does 
not depend strongly on local density. 

There are various scenarios for why the field and 
cluster fractions mi ght be the same. As described in 
IMartini et al.l (|2009f) . it is possible that blue, late- type 
galaxies in clusters are still in-falling and have not yet 
been stripped of their gas by interactions with the clus- 
ter potential. As a result, these galaxies might retain 
enough cold gas to fuel an AGN, as is the case in the 
field. In another scenario, red, early- type galaxies may 
keep their cold gas (despite interactions with cluster po- 
tential) and thus continue to host an active SMBH. A 
third possibility invokes evolution of the AGN fraction 
in both clusters and in the field — here, AGN activity 
was more dependent on environment in the past, e.g., 
closer to the peak in the AGN number density a,t z ^ 2, 
but at the present day galaxies in both environments are 
relatively quiescent and the differences between the field 
and clusters is not discernible. This last scenario can be 
tested by probing the AGN fraction in both clusters and 
the field at high redshifts (and possibly in regimes where 
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low-luminosity AGN dominate). 

As discussed in and ^16. 11 there is some evidence 
that the AGN fraction evolves as a function of redshift. 
Though we find the AGN fraction in clusters and in the 
field in close agreement at low-redshift, they may begin 
to differ as we probe to hi gher r e dshift regimes. Em- 
ployin g a method similar to iPagd (|2Q01[ ). iMartini et al.l 
()2009() estimate the field AGN fraction based on galaxy 
luminosity functions and compare it to their Chandra 
cluster fraction (see their Fig. 3). At z ~ 0.15, they find 
the two in agreement, but at z ~ 0.8 the field fraction is 
nearly a factor of 5 higher than the cluster fraction. This 
gives some credence to a scenario wherein the cluster and 
field AGN fractions were more disparate at ea rlier times . 

Other authors (jGeorgakakis et all l2008al : 
[Silverman et al.l l2009al ) have investigated the AGN 
fraction in groups of galaxies (or as a function of local 
over-densities). Using DEEP2 Chandra observations 
from AEGIS for 0.7 < z < 0.9 and Mb < -20, 
iGeorgakakis et al.l (|2008al) find group and field fractions 
of 4.7 ± 1.6% and 4.5 ± 1.0%, respectively. This group 
fraction is ab o ut 5 t imes higher than that reported by 
IMartini et al.l (|2009[ ) in clusters at similar redshifts. 
This discrepancy may arise from intrinsic differences in 
the group and cluster environments at redshifts near 
z ^ I, but more likely it originates from differences in 
the X-ray and optical se lection criteria of the surveys; 
i.e. IGeorgaka kis et aTl (l2008af) include l ower- X -ray- 
luminosity sources than do iMartinf "eTall (pool ) and 
their host galaxy magnitude limits d i ffer. U sing only 
the AEGIS data, IGeorgakakis et al.l ()2008aD report a 
higher AGN fraction in groups than in the field at 
z ^ 1, but attribute this to their finding that AGN are 
preferentially found in red, luminous galaxies, which 
are themselves more likely to reside in higher density 
environments. They note that an exception to this trend 
may exist for the most X-ray-luminous AGN, which 
a ppear in higher numbers in the field. 

iSilverman et al.l (|2009aD also study the X-ray AGN 
fraction for low-mass host galaxies (10.4 < log M^, < 11) 
in the field and in groups, and find no statistical dif- 
ference between the fractions in these environments. 
They do report a trend for log M^, > 11, wherein 
high X-ray luminosity A GN are less common in more 
dense environs, s ee als o iKauffmann et al.l ()2004[ ) and 
iConstantin et al." (200§). Fo r 0.58 < z < 1.05 and log 
> 10.4, Silverman et al.l (|2009aD find a field fraction 
of ^ 3.7 ± 0.5% and a grou p fraction of ^ 4zfc 1 2%, c on- 
sistent with the findings of IGeorgakakis et al.l (|2008aD . 

We look to our highest redshift bi ns to compare our 
field X-ray AGN fractions to t hose o f IGeorgakakis et"al] 
(p008a) and ISilverman et al.l (|2009al ). For 0.55 < z < 
0.7, our only complete mass bin is 11.6 < log Ai* < 13, 
for which we calculate an upper limit on the AGN frac- 
tion of 25.00 ± 20.80%. Note, however, that we base this 
estimate of the fraction on only two X-ray-observed host 
galaxies and zero X-ray detections and thus do not of- 
fer a sensitive probe of the field fraction (particularly 
within this redshift range). For 0.4 < z < 0.55, we 
have two complete mass bins, 11 < log < 11.6 and 
11.6 < log M^, < 13, in which we measure AGN fractions 
of 5.33±1.70% and 3.70±3.13%. Though these fractions 
are broadly consistent with the AEGIS and zCOSMOS 
studies, they sample higher mass hosts at lower redshifts 



— our median redshifts are {z) = 0.44 and (z) = 0.45 
for the low- and high-mass bins, respectively. To include 
lower-mass hosts (e.g., log iW^* > 10.4), we must set our 
redshift limit to z = 0.4 to achieve a complete sample, 
and questions about evolution begin to impede compar- 
isons with the deep field results. 

We point out these higher-redshift results to illustrate 
that small number statistics continue to stymie measure- 
ments of the X-ray AGN fraction in all environments 
(field, group, and cluster), particularly at the faint end 
of the galaxy population and at high- redshift. 

6.3. AGN Host Calaxy Properties 

To investigate the properties of our host galaxies, we 
show the AGN fraction in bins of bestz vs. restframe 
(u— r)o and {U—V)o color (Fig. [121 left and right panels, 
respectively). The (m — r)o blue, green, and red color bins 
are (u-r)o' < 1.8, 1.8 < (u-r)o < 2.6, and (u-r)o > 2.6. 
The green bin straddles the (u — r) p = 2.2 divide betwee n 
blue and red galaxies defined by iStrateva et al.l (|2001| ) . 
The (U - V)o color bins are 0.8 < {U - V)o < 1.4, 
1.4 < {U - V)o < 1.8, 1.8 < {U - V)o < 2.1, and 
2.1 <{U - y)n < 2.4, selected to match the zCOSMOS 
color bins ^Silverman et al.l l2009b. their Figure 12b). We 
apply a mass cut of log Mi, > 10.6 (black solid line in 
the right panel of Fig. [9]), again to match the zCOSMOS 
study of host galaxy colors, and include only the redshift 
ranges associated with our samples l*-3*, as these are 
complete for this mass cut. For the (u — r)o comparison 
we employ the X-ray limit used throughout this work 
{Lx = lO''^ ergs-i); for the zCOSMOS {U - V)o com- 
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ergs ^ as in lSilverman et all 



parison we use Lx 
(|2009bD . 

In the {u — r)o bins, the fraction of AGN in the 
blue sequence is larger than that in cither the red se- 
quence or the green valley, for all three redshift intervals 
(Fig. [ini left panel, colored symbols). This trend to- 
ward a larger AGN fraction in the blue sequence persists 
when the red point sources, which might bias the red 
sequence and/or the green valley toward a lower AGN 
fraction, are removed (Fig. \E\ green points; see also 
i 35.5p . "Upper limit" AGN fractions, calculated with- 
out these red point sources, are shown in Figs. [T3] and 
[TBI (filled, black, downward- facing triangles). This result 
is in contrast to multiple studies, which have found that 
AGN lie preferentially in the valley between the red and 
blue galaxy sequ ences, e.g., Georgakakis et al, ( 2 008a. b): 



Silverman et al 



Silverman et al 



2008); Sch awinski et al.l (|20ia ) : though 
2008 ) associate the green valley peak 
with large-scale structures at particular redshifts. Our 
finding agrees, however, with more recent studies of the 
zCOSMOS fields (Silverman et al. 2009a. 1?); these au- 
thors argue that the peak reported in the transition re- 
gion between blue and red galaxies is an artifact that can 
be eliminated by applying an appropriate mass cut, as 
we have done here. 

As discussed in the AGN central engine may con- 
tribute considerable light to the SDSS SED, making the 
object appear artificially blue. AGN contamination of 
the host galaxy light would result in an overestimate of 
the number of galaxies along the blue sequence. We do 
not expect this effect to be large, however, since the ma- 
jority of the X-ray sources in our three lowest redshift 
bins {z < 0.4) have Lx < 10^" ergs"! (Fig. M left 
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42, solid black line in the right panel of Fig. [gj for three low-redshift samples, binned 

color and the 
to compare to 



)o color (see ^for discussion). (Right) Similar to the left panel, but the data are binned instead by {U — V_)o 

scale (black dashed line in the right panel of Fig. [6j 



L'y = 42.5 and plotted on a \ 



zCO SMO S results (open red diamonds). Upper limits on the fractions are plotted as solid, downward- facing black triangles in both panels 
(see i|5.5ll . The trend toward lower AGN fractions at redder colors is clearly visible for both color regimes. 
(A color version of this figure is available in the online journal.) 

panel); hence, accretion onto the central SMBH is prob- 
ably not powerful en ough to sig nificantly contaminate 
the host galaxy light (|Prieto et a l. 2010). A population 
of starburst galaxies, whose X-ray luminosities can be as 
high as Lx ~ 10''^ ergs~^, might also contaminate our 
blue sequence, but such X-ray-bright star-forn iing galax- 
ies ar e more likely at higher redshift {z > 0.4: lYan et al.l 

IMoh . 

To further explore the properties of our sample as 
a function of host galaxy properties, we compare our 
primarily photometric Full galaxy sample to our spec- 
troscopic sub-sample. Figure [17] shows {u — r)o color 
vs. absolute r-band magnitude (M^; left panels), and 
{u — r)o color histograms (right panels) for each of our 
five volume-limited samples. In the left-hand panels we 
overplot spectroscopic types for all galaxies that have 
secure line identifications in the MPA/JHU DR7 cata- 
log, as well as those identified via ChaMP spectroscopic 
follow-up (filled symbols indicate X-ray detections, open 
symbols mark X-ray non-detections); blue squares are 
narrow emission line galaxies (NELG); red triangles show 
absorption line galaxies (ALG); black circles mark broad- 
line AGN. We also indicate X-ray-detected sources, with 
and without spectroscopic classifications, as open green 
diamonds. In the right-hand panels, we show (u — r)o 
histograms for the Full galaxy sample (solid black line), 
the X-ray-detected sub-sample (green dotted line), as 
well as for the MPA/JHU NELG (blue dashed hue) and 
ALG (red dot-dashed line) samples. The histograms are 
each normalized to the total number of objects within 
the (sub-)sample, shown in parentheses in the relevant 
figure. 

We see a clear division between the NELG and ALG 



samples in Fig. [T71 with NELGs falling along the 
blue sequence and ALGs along the red sequence. We 
do not see an over-density of X-ray sources in the 
green valley between the red and blu e sequences, as has 
somet im es been reported elsewhere (iGeorgakakis et al.l 
[2008allbl : [Sih^erman et al.l [20081 : ISchawinski et al.ll2010D . 
In the color histograms, the Full galaxy sample demon- 
strates a clear bi-modality at low redshift, but this 
trend becomes washed out (or perhaps shifted) at higher- 
redshift. This may result from evolution in the galaxy 
population, from smearing due to photo-z uncertainties, 
or from stellar interlopers. As an illustration of this last 
possibility, we note that for Sample 1 (Fig. [T71 first 
panel) the peak of the red sequence in the Full galaxy 
sample (solid black line) is redder than the spectroscopic 
ALG sample (red dot-dashed line) — this offset may 
be due to a population of M-dwarf stellar contaminants 
amongst the objects with photometric redshifts. Here 
again, contamination of the host color by a luminous 
central engine may force galaxies out of the green or red 
portions of the diagram and onto the blue sequence (see 
§5.5|) . blurring the expected bi-modality in the underly- 
ing galaxy population. However, since such contamina- 
tion is unlikely in the majority of the non- X-ray-detected 
(Full galaxy) sample, we do not expect the effect to be 
large enough to explain the trends we see here. 

We also compare the {U — V)o results (Fig. [TH 
right panel) to the field fraction determined in the 10k 
catalog of the zCOSMOS spec troscopic redshift survey 
([Silverman et a l. 2009a^lJ). The [Silverman et al.l ([2009b[ ) 
fractions for these four {U — V)o bins (estimated from 
their Fig. 12b) are approximately 3.5 ±0.9%, 3.9 ±0.6%, 
1.2 ± 0.3%, and 0.25 ± 0.35% (Fig. [H right panel, red 
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Fig. 17. — {Left sub-panels) Restframe (u — r) as a function of r-band absolute magnitude for the Full galaxy sample (black dots), the 
X-ray-detected sample (open green diamonds), and objects with spectra: ChaMP/SDSS NELG (open blue squares), ChaMP/SDSS X-ray- 
detected NELG (filled blue squares), ChaMP/SDSS ALG (open red triangles), ChaMP/SDSS X-ray-detected ALG (filled red triangles), 
SDSS BLAGN (open black circles), and ChaMP BLAGN (filled black circles). (Right sub-panels) Histogram of (n - ro) for the Full 
photometric sample (black solid line), the X-ray-detected sample (green dotted line), the SDSS spectroscopic NELG sample (dashed blue 
line), and the SDSS s pectroscopic ALG sample (dot-dashed red line). The (n — r)o = 2.2 division between the red and blue sequences from 
IStrateva et al.l 120011) is shown in both panels (dashed black line). 
(A color version of this figure is available in the online journal.) 



open diamonds). We compare these to our {U — V)o 
fractions (taking a weighted averaged over our three red- 
shift bins) of 0.67 ± 0.12%, 0.26 ± 0.10%, 0.14 ± 0.08%, 
and 0.40 ± 0.19%. As was the case for the global field 
fraction (gUI]), our ChaMP/SDSS fractions faU signifi- 
cantly lower than the zCOSMOS fractions, except in the 
reddest color bin. The zCOSMOS sample is almost cer- 
tainly weighted toward higher redshifts (their redshift 
range is 0.1 < z < 1.02, vs. ours for this compari- 
son, 0.0025 < z < 0.4). Thus, we might rather com- 
pare only to our AGN fraction in the highest redshift 
bin (0.275 < z < 0.4), i.e. 0.54 ± 0.19%, 0.93 ± 0.37%, 
0.51 ±0.33%, and 0.36 ±0.3%. Even here, our field AGN 
fractions are systematically lower than those found in the 
zCOSMOS survey; they disagree by approximately 3 a 



in the blue and green sequences, by roughly 2 cr in the 
red sequence, and agree within the errors only in the red- 
dest color bin. Despite these differences, both our AGN 
fractions and the zCOSM OS fractions are la r gest fo r host 
galaxies with blue colors. [Silverman et al.l (l2009bO pro- 
pose that this trend might establish a connection between 
accretion onto SMBH and star formation out to z '--^ 1. 
Our findings lend additional support to this hypothesis. 

It is plausible that the difference between the field 
AGN fractions from ChaMP/SDSS and zCOSMOS arises 
from evolution in the AGN fraction. If we assume a mean 
redshift for the zCOSMOS sample of z ^ 0.6, and a mean 
redshift of z ^ 0.3 for our sample 3, we would expect PLE 
(1 + zf evolution in the AGN XLF to yield a zCOS- 
MOS fraction that is a factor of ^ 2 higher than the 
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ChaMP/SDSS fraction. Such a correction would bring 
our results closer to the zCOSMOS blue and green se- 
quence fractions, and achieve agreement in the red se- 
quence. Note that the fractions in the reddest bin would, 
however, begin to disagree. Evolution of the underlying 
host galaxies colors may explain these differences, e.g., 
AGN activity might occur most often in blue galaxies at 
higher redshifts {z > 0.5), but increasingly shift toward 
red galaxies in the local Univ erse (see also discussions in 
[Silverman et"anilM 120090 ). 

7. CONCLUSIONS 

We present a study of the X-ray-active fraction of field 
galaxies in an extensive catalog of more than 100,000 
galaxies with X-ray and optical coverage from Chandra 
and SDSS (^1,600 are X-ray-detected). We combine 
ChaMP/SDSS spectroscopic and photometric redshifts 
with X-ray and optical fluxes (or flux limits) to assign ab- 
solute magnitudes. X-ray luminosities (or limits), masses, 
and colors to each galaxy in the sample. With these 
data we explore the AGN fraction {Lx (0.5 — 8keV) > 
j^q42 ergs~^) in five independent samples, complete in 
redshift and i-band absolute magnitude: 

• Fagn.i = 0.16 ± 0.06% {z < 0.125, -18 > K'h > -20) 

• Fagn,2 = 0.50 ± 0.11% {z < 0.275, -20 > Mi > -21) 

• Fagn,3 = 1.27 ±0.18% {z < 0.400, -21 > Mi > -22) 

• Fagn,4 = 2.85 ± 0.39% {z < 0.550, -22 > Mi > -23) 

• Fagn,5 = 3.80 ± 0.92% {z < 0.700, Mi < -23) 

The low-redshift bin edge for these fractions is z = 
0.0025, selected to avoid galactic objects and photomet- 
ric redshift artifacts. 

Our robust analysis is enabled by ChaMP's compre- 
hensive sensitivity maps for ACIS imaging, which al- 
low recognition of imaged-but-undetected objects, counts 
limits for 50% and 90% detection completeness, and cor- 
responding flux upper limits at any sky position. The 
AGN (or X-ray-active) fraction for each volume-limited 
sample, as well as the asymmetric error bars, are calcu- 
lated via the /3 distribution — an ideal statistical tool 
for evaluating fractions, particularly when the numera- 
tor and denominator are not drawn from independent 
samples. 

We find excellent agreement between our 
ChaMP/SDSS field AGN fraction and the Chandra 
cluster AGN fraction, for samples restricted to sim- 
ilar redshift and absolute magnitude ranges. For 
ChaMP/SDSS field galaxies with 0.05 < z < 0.31 and 
absolute i?-band magnitud e more luminous tha n -20, 
Fagn = 1.19 ± 0.11%. The lMartini eFall (12001 X-ray 
cluster fraction for this redshift and Mn is 1.0%. We 
find good agreement between our field fraction and the 
cluster fraction for a variety of other absolute i?-band 
magnitude and X-ray luminosity limits as well. Our 
results are also broadly consistent with measures of 
the field AGN fraction in narrow, deep fields, though 
differences in the optical selection criteria, redshift 
coverage, and possible cosmic variance between fields 
introduce larger uncertainties in these comparisons. 

In our analysis of the AGN fraction, as well as in our 
comparison to deep field studies, we find evidence that 
the AGN fraction evolves with redshift. Our data are 
consistent with either (1 -|- z)^ or (1 -I- z)"^ evolution, i.e. 



the two most common fits to the X-ray AGN luminos- 
ity function, but likely also depend on evolution in the 
host galaxy population. These findings are tantalizing, 
but poorly constrained and require high-quality wide- 
field data at redshifts out to (and above) z ~ 1 for ver- 
ification. We test the impact of host galaxy color (for 
stellar masses log M^, > 10.6) and find that the AGN 
fraction is largest for hosts with blue colors, though this 
trend too may depend on redshift. Since galaxy and 
AGN evolution are both affected by environment, we are 
embarking on a clustering study of AGN activity vs. en- 
vironment to quantify these properties. A comparison 
between cluster, group, and field fractions at redshifts 
approaching z — 2 would be particularly compelling. 
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Facilities: CXO, Sloan, FLW0:1.5m (FAST spectro- graph), Magellan:Baade (LDSS2 imaging spectrograph), 

MagehaniClay (IMACS), WIYN (Hydra) 

APPENDIX 

In we define the fraction of X-ray-active galaxies. We apply both optical and X-ray completeness criteria to 
insure an accurate measure of the fraction. In this Appendix, we describe the objects that appear in our volume-limited 
optical samples, but do not pass the X-ray cuts. These "dropped" objects can be described by: 

Zmax A'^bright OO 

H ^G{LxMm,M„bestz) , (1) 

Ndet,drop = ^^X{Lx,Lx,liTn,Mi,bestz) , (2) 

where Num^drop quantifies the number of X-ray-observed objects "dropped" from the denominator, and Ndet,drop 
describes the number of X-ray-detected sources eliminated from the numerator. 

In Fig. [THl we show the behavior of Ndet, Num, Ndet,drop, and Num^drop as a function of limiting X-ray luminosity 
for sample 4. Comparing this to the top, right panel of Fig. [TOl it is clear that the high X-ray-active fraction below 
log L'-^ ^ 41.4 can be explained by a relatively small number of both X-ray-detected and X-ray-observed galaxies. At 
higher limiting X-ray luminosities, the number of X-ray detections in the numerator increases (until log L'^ 43), but 
the number of galaxies that might have been detected at this threshold X-ray luminosity (the denominator) increases 
even more rapidly, leading to an overall decline in the X-ray-active fraction. At the highest luminosities the number 
of X-ray-detected sources falls off sharply while the number of X-ray-observed galaxies continues to climb; the fraction 
drops to zero, accordingly. The error intervals shown in Fig. [10] are intuitive insofar as they are largest where the 
number of objects used to determine the fraction is small, and smallest in the opposite extreme. 
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Fig. 18. — The numbers of X-ray-detected (A^det) and X-ray-observed {Nn^) galaxies included in the X-ray-active fraction at each X-ray 
threshold for sample 4 are shown as black dashed and solid lines, respectively. The numbers of objects that are "dropped" at each limit 
are also shown (gray dot-dashed line for X-ray-detections \Ndet.dro v\ i g ray dotted line for X-ray-observed objects [Nu^^drop])- This plot 
makes explicit the shape of the X-ray active fraction curve in fig. 1101 i.e. while the number of X-ray-detections included in the sample 
ranges from a few to several tens of objects between log Lx of 40.5 and 44, the number of X-ray-observed galaxies increases from about 
20 to more than 20,000 in the same range. In this optically-complete sample the majority of the brightest X-ray sources make it into the 
fraction, while many X-ray-low-luminosity sources are eliminated. 
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